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ABSTRACT 

We present a Chandra study of mass profiles in 7 elliptical galaxies, of which 3 have galaxy-scale and 
4 group-scale halos, demarcated at 10 13 M©. These represent the best available data for nearby objects 
with comparable X-ray luminosities. We measure ~flat mass-to-light (M/L) profiles within an optical 
half-light radius (R c ff), rising by an order of magnitude at ^10R o ff, which confirms the presence of 
dark matter (DM). The data indicate hydrostatic equilibrium, which is also supported by agreement 
with studies of stellar kinematics in elliptical galaxies. The data are well-fitted by a model comprising 
an NFW DM profile and a baryonic component following the optical light. The distribution of DM 
halo concentration parameters (c) versus M V j r agrees with ACDM predictions and our observations 
of bright groups. Concentrations are slightly higher than expected, which is most likely a selection 
effect. Omitting the stellar mass drastically increases c, possibly explaining large concentrations 
found by some past observers. The stellar M/Lk agree with population synthesis models, assuming 
a Kroupa IMF. Allowing adiabatic compression (AC) of the DM halo by baryons made M/L more 
discrepant, casting some doubt on AC. Our best-fitting models imply total baryon fractions ^0.04- 
0.09, consistent with models of galaxy formation incorporating strong feedback. The groups exhibit 
positive temperature gradients, consistent with the "Universal" profiles found in other groups and 
clusters, whereas the galaxies have negative gradients, suggesting a change in the evolutionary history 
of the systems around M V j r ~ 10 13 M Q . 

Subject headings: Xrays: galaxies — galaxies: elliptical and lenticular, cD- 
ISM— dark matter 



galaxies: halos — galaxies: 



1. INTRODUCTION 

The nature and distribution of dark matter (DM) in 
the Universe is one of the fundamental problems facing 
modern physics. Cold DM lies at the heart of our current 
(ACDM) cosmological paradigm, which predicts substan- 
tial DM halos for objects at all mass-scales from galax- 
ies to clusters. Although ACDM has been re markably 
successful at explaining large-sca le features (e.g. Spergel 
|st al. [2003 ; Pcrlmutter ct aL 199S ) , observations of galax- 
ies have been more problematical for the theory. Dissi- 



mological parame ters, including a* a nd the dark energy 
equation of state ( Kuhlcn ct al. 2005| ), making an obser- 
vational test of this relation potentially a very powerful 
tool for cosmology. 

The mass profiles of galaxies also may provide valu- 
able clues as to the way in which galaxies form in DM 
halos. In particular, as baryons cool and collapse into 
stars, the associated increase in the central mass den- 
sity should in turn modify the shape of the DM halo. 
This process is typically modelled assuming adia batic 



patiq nleaa dark matter aimulationa find that dark matter 



contraction (AC) of the DM particle o rbits (e.g. Blu- 
hai^ ^ mcnthal ct al.|[l98(| |Gnedin etalj jooj). If the gaW 



sity profile (Navarro et al. 1997, hereafter NFW) over a 



wide range of Virial masses (M v ; r ) (e.g. [Bullock et al 
[2001 ). Low mass halos tend to form first in hierarchi- 
cal cosmologies and are consequently more tightly con- 
centrated than their later forming, high mass counter- 
parts. This tendency produces a predicted correlation 
between the DM halo concentration parameter (c, which 
is ratio between Virial radius, R v ir, and t he characteris- 
tic sc ale of the density profile) and M v ; r ( Navarro et al. 
1997). However, since mass and formation epoch are not 
perfectly correlate d, we expect a significant scatter at 



halo subsequently evolves by major mergers, simulations 
are unclear as to whet her these features would persist 
(e.g. Gncdin ct al. 2004) or whether the merging process 
may destroy this imprint of star formation, or even mix 
the DM and baryons sufficiently to prod uce a total grav- 
itating mass profile more akin to NFW (Loeb & Peebles 



2003; El-Zant et al. 2004) 



fixed Virial ma ss ( [lingj P000| ; [Bullock et al.||2001| ; |Wech 



lcr ct al. 2002). The tight link between halo formation 



epoch and concentration implies that the precise relation 
between c and M v j r is sensitive to the underlying Cos- 

nRpartmpnt nf Prysirs anH Astrnnnmy University nf fla.lifnr- 



Observational tests of the predictions of ACDM have 
proven controversial. In clusters of galaxies there is over- 
whelming evidence for DM, and an increasing body of 
work verifying the predictions of ACDM. In particular re- 
cent, high-quality Chandra and XMM observations have 
revealed mass profiles in remarkable agreement with the 
Universal profile f rom deep in the c o re to a large frac 



tion of Ryj r (e.g. |Lcwis ct al. 2003 ; Zappacosta et al 



2006; Vikhlinin ct al. 200£), and a distrib ution of c 



sus M v j r in good agreement with ACDM (Pointecouteau 
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nia at Irvine, 4129 Frederick Reines Hall, Irvine, CA 92697-4575 Ct al. 2005). In galaxies, however, the picture is much less 



clear. Rotation curve analysis of low surface brightness 
(LSB) disk galaxies has suggested si gnificantly less cusp y 
density profiles than expected (e.g. Swatcrs et al. 2000] ) . 
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Although this discrepancy led to a serious dis cussion of 
mod ifications to the standard paradigm (e.g . Hogan & 
Dale mton |2000|; ISpcrgcl fe Steiuhardtl |200C|; gentner fc| 



Bulk[ck| |2002| ; |Kaplinghat| |2005| ; |Ccmbranos et al.| p005[ ) 

recent results, taking account of observational bias and 
the 3-dimensional geometry of the D M halos. have done 



relatively rare instances of strong lensing are required to 
study DM in individual systems. Nonetheless it has been 
possible in a few cases to decompose the mass into stel- 
lar and DM components, albeit with strong assumptions 



much to reso l ve th e discrepancy (e.g. S waters et al. 2003 



|5imon et al]|2005| ). However, some significant discrep- 



or additional observational con straints (e.g. Rusin et al 
20031 ; [Treu 
X^ray observations 



Koopmans pb04| ) 

is of the hot gas in early-type galax- 



ancies remain, not least of which is that the DM halos 
of th ese galaxies appear less concentrated th an expected 
(e.g. Gonzalez et al. 2000 ; Kassin et al. 2006| ). A possible 



explanation is that LSB gala xies are preferentially found 
in lo w-concentration halos (Bullock et al] 2001; Bailin 
et al. 2005 ; Wechsler et al. 2005 ) , making additional con- 
straints at the galaxy scale extremely important. 

In many respects, kinematical mass measurements are 
far more challenging for early- type than spiral galaxies. 
As essentially pressure-supported systems little is known 
a priori about the velocity anisotropy tensor of the stars 
in elliptical galaxies, which is problematical for the deter- 
mination of the mass from stellar motions. Nonetheless, 
stellar kinematical measurements have widely been used 
as a means to measure the gravitating matter within 



ies provide a complementary means to infer the mass 
profiles via techniques similar to those used in study- 
ing clusters. Since the X-ray emission from early- type 
galaxies is typically not very bright, prior to the advent 
of Chandra and XMM this was limited by the relatively 
sparse information on the radial temperature and den- 
sity profiles of the hot gas which could be determined 
by prior generations of satellites. Notwithstanding this 
limitation, large M/L ratios (consistent with substantial 
DM) were inferred for a number of X-ray bright galax- 
ies, albeit with strong assumptions concerning the tem- 
perature and density profile s (e.g. Forman et al.| 1985 



Loewenstein fc White 1999 ). Using a novel technique 
which relied, instead , on th e ellipticitv of the X-ray halo, 
Buote fc Canizares ( 1994 ) were able robustly to detect 



the optical half-light radius (R c ff) of elliptical galax- 

. ell' 



DM in the isolated el l iptical NGC 720 (s ee also [Buote 



ies ( e.g. Bmncy et al. 199C; van der Mareil 119911 IGer 
hard et aT. 2001 ). These studies tend to find relatively 
flat mass-to- light (M/L) ratios within R ff , implying that 
most of the matter within this radius is barvonic. Con- 



Canizares||1996| , |1998j ; |Buote et al.||2002| ). Detailed mea- 
surements oFEhe ramalnMss - a[istribution were, however, 
largely restricted to a few massive systems, which may be 
entwined with a group halo (e.g. Irwin fc Sarazin 199€ 
|BrighciitT 



sider ition of the tilt in the fundamental plane can also Mathews \ 1997 ) were able to decompose the mass profiles 



20031). In 



lead to t he same conclusion flBorricllo et al, t 
contrast, Padmanabhan et al.] ~( 2*0'04| ) pointed out that 
dynamical M/L ratios within R ff are much larger than 
predicted by realistic stellar population synthesis models 
for stars alone, allowing > 50% of the mass within R e g 
to be dark. 

Attempts to extend kinematical studies of elliptical 
galaxies to larger radii, where DM should be dominant, 
have p roven controversial. In particular Romanowsky 
et al] ( 2003| ) argued against the existence of DM in a 
small sample of elliptical galaxies, based on planetary 
nebulae dynamics within ^5R c ff . We note that this sam- 
ple was heavily biased towards very X-ray faint objects, 
which might hint at low-mass halos since they have no t 
held onto their hot gas. In any case Dekel et al. ( 2005| ) 
pointed out that their conclusions were very sensitive 
to the uncertainty in the velocity anisotropy tensor, for 
plausible values of which the data were consistent with 
substantial DM halos. In fact globular cluster dynamics 
in one of these system s, NGC 3379 , does imply a signif- 



Mathews 1997). Nevertheless Brighcnti 



icant amount of DM (Pierce et al. 2006; Bergond et al 
p006| ). As more kinematical studies of early- type galaxies 
at large radii are appearing, it is becoming clear that at 
least some elliptical gal a xies host considerab le DM halos 
(e.g. 3tatler et al. 1999; Romanowsky 2005). There per- 
sist some questions, however, as to the extent to which 
all galaxie s have DM halos c onsist ent with ACDM. In 
particular Napolitano et al. (2005) argued that a sub- 
stan tial number of early-type galaxy halos appear less 



of two systems, NGC 4472 and NGC 4 649, into stellar 
and DM components. |Sato et al ( 200C| ) investigated the 
M v i r -c relation using ASCA for a sample of objects rang- 
ing from massive clusters to ^3 elliptical galaxies. The 
limited spatial resolution of ASCA necessitated some as- 
sumptions about the density profiles and, crucially, the 
authors neglected any stellar mass component in their 
fits. This omission may explain the very steep M v ; r -c 
relation (with 0200^ 30 for the g alaxies) found by thes e 
authors, in conflict with ACDM flMamon fc Lokas||2005| ). 

Although mass profiles of early-type galaxies are be- 
ginning to appear which exploit the improved sensitivity 
and resolution of Chandra and XMM, many of the most 
interesting constraints on DM are still restricted to mas- 
sive syste ms, which may b e at th e centres of groups. For 
example, Fukazawa et al.| ( [2006 ) reported Chandra and 
XMM M/L profiles for ^50 galaxies and groups, confirm- 
ing ^flat profiles within R ff which rise at larger radii. 
However, the constraints at large radii were dominated 
by the massive (group-scale) objects so the implications 
for the DM content of normal galaxies are unclear. Fur- 
thermore, the authors included a substantial number of 
highly disturbed systems, in which hydrostatic equilib- 
rium may be questioned, and failed to account for the 
unresolved sources which dominate the emission in the 
lowest-Lx objects in their sample 4 . Recently, however, 
detailed Chandra and XMM mass profiles have begun 
to appear for isolated early-type galaxies, also confirm 



conci :ntratcd than expected 



ing the presence of massive DM halos ( e.g. p' Sullivan 
Ponman|[20"ol |Khosroshahi et al.||2004[ ). 



Gravitational lensing provides further evidence that, at This paper is part of a series (see also |Gastaldello et al 



Kochanek 

1 rp - ' 


1995); Fischer et al. 


200C ; 


Rusin 



vides useful mass constraints in a statistical sense, the 



4 Although the authors account for unresolved sources when 
measuring the gas temperature, they do not account for it when 
computing the gas density, where its effect is more pronounced 
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[2006 ; Zappacosta et al. 2006 ; Buotc ct al. 2006a| ,|b|) us- 
ing high-quality Chandra and XMM data to investigate 
the mass profiles of galaxies, groups and clusters. This 
provides an unprecedented opportunity to place defini- 
tive constraints upon the M V j r -c relation over ~2 orders 
of magnitude in M v ; r . In this paper, we focus on the 
temperature, density and mass profiles of seven galax- 
ies and poor groups chosen from the Chandra archive. 
In order to compare to theory we perform spherically- 
averaged analysis, leaving a discussion of the ellipticities 
of the X-ray halos to a future paper. In § || we discuss 
the target selection. The data-reduction is described in 
§ H and the X-ray morphology is addressed in § 0. We 
discuss the spectral analysis in § ^, the mass analysis in 
§ |^, the systematic uncertainties m our analysis in § 
and reach our conclusions in § ||. The three systems for 
which we find M V j r < 10 13 M Q are optically isolated and 
so we refer to them as "galaxies" , and the other systems 
in o ur sa mple as groups. We discuss this in more detail 



undergone a major merger. For the purposes of this pa- 
per, however, we do not require statistical completeness, 
and we will discuss how to take account of these selection 



3.3 



In this paper, all error-bars quoted represent 
90% confidence limits, unless otherwise stated, and we 
computed Virial quantities assuming a "critical overden- 
sity" factor for the DM halos of phaio/ Pent = 103 (where 
Phaio is the mean density of a sphere of mass M v ; r and 
radius R V ir)- 

2. TARGET SELECTION 

We chose, for this initial study, to focus on objects 
observed with Chandra. Chandra data are particularly 
valuable for the study of galaxies since the unprecedented 
spatial resolution makes it possible to resolve the temper- 
ature and density profiles deep into the galaxy core, al- 
lowing us to disentangle the stellar and dark matter, and 
resolve them into discrete components. We initially chose 
a set of potential target systems from det e ction s listed 
in the X-ray catalogue of O'Sullivan et al. ( 2001 ) which 
have non-grating ACIS data in the Chandra archive. To 
eliminate bright groups and cluster cDs in the sample, 
we excluded galaxies with Lx^ 10 43 erg s _1 . In order to 
perform the required spatially-resolved spectroscopy, we 
required at least ~5000 hot gas photons. The potential 
targets were processed and the 0.1-10.0 keV image ex- 
amined for evidence of large-scale disturbances (§ 0). We 
included some systems with low-amplitude asymmetries 
which should not strongly disturb hy dro static equilib- 
rium (we discuss this in more detail in § |7.4| ) . Preliminary 
analysis was conducted to estimate the Virial mass of the 
object (§ [|). Since we aimed to focus on lower-mass ob- 
jects, systems for which a fit using a simple NFW profile 
yielded M v ; r > 10 13 M Q were discounted. Mas sive objects 
of this type are the focus of another study ( Gastaldcllo 



|st al 



2006 ). The most promising candidates for study 
found via this method were chosen for detailed analysis. 
The properties of the 7 objects in our sample and the 
Chandra exposures are shown in Table |l|. 

Our selection criteria naturally bias the sample towards 
X-ray bright galaxies. One might expect that galaxies 
sitting in deep potential wells are more likely to retain 
hot gas than those with little dark matter, and so our 
results may be biased somewhat towards those galaxies 
with substantial dark halos fin contrast to the opposite 
bias in the analysis of Romanowsky et al. 2003 ). As we 
are selecting objects which are not heavily disturbed, we 
are also biased towards galaxies which have not recently 



effects in Buote et alj (2006a 



3. DATA REDUCTION 

For data reduction, we used the CIAO 3.2.2 and Hea- 
soft 5.3 software suites, in conjunction with Chandra cal- 
ibration database (Caldb) version 3.1.0. Spectral- fitting 
was conducted with Xspec 11.3.1w. In order to ensure the 
most up-to-date calibration, all data were reprocessed 
from the "level 1" events files, following the standard 
Chandra data-reduction threads 5 . We applied the stan- 
dard correction to take account of the time-dependent 
gain-drift as implemented in the standard CIAO tools. 
To identify periods of enhanced background ("flaring"), 
which seriously degrades the signal-t q- noise fS/N) and 
complicates background subtraction ( Markevitch 2002 ) 
we accumulated background lightcurves for each expo- 
sure from low surface-brightness regions of the active 
chips. We excluded obvious diffuse emission and data 
in the vicinity of any detected point-sources (see below) . 
Periods of flaring were identified by eye and excised. 
Small amounts of residual flaring not removed by this 
procedure can be important in low surface-brightness re- 
gions at large radii, but this was take n into account in our 
treatment of the background (§ 3.1). The final exposure 
times are listed in Table |l[ 

Point source d etection was perform ed using the CIAO 
tool wavdetect (Freeman ct al. [2002 ). Point sources were 
identified in full-resolution images of the ACIS focal- 
plane, containing all active chips (except the S4 chip, 
which suffers from serious "streaking" , which can lead to 
false detections). To maximise the likelihood of identify- 
ing sources with peculiarly hard or soft spectra, images 
were created in three energy bands, 0.1-10.0 keV, 0.1- 
3.0 keV and 3.0-10.0 keV. Sources were detected sep- 
arately in each image. In order to minimize spurious 
detections at node or chip boundaries we supplied the 
detection algorithm with exposure-maps generated at en- 
ergies 1.7 keV, 1.0 keV and 7 keV respectively (although 
the precise energies chosen made little difference to the 
results). The detection algorithm searched for structure 
over pixel-scales of 1, 2, 4, 8 and 16 pixels, and the detec- 
tion threshold was set to ensure ~0.1 spurious detections 
per image. The source-lists obtained within each energy- 
band were combined and duplicated sources removed, 
and the final list was checked by visual inspection of the 
images. The data in the vicinity of any detected point 
source were removed so as not to contaminate t he diff use 
emi ssion. As discussed in Humphrey fc Buote ( 2004 , see 
also Kim k, Fabbiano 2004 ) a significant fraction of faint 
X-ray binary sources will not have been detected by this 
procedure, and so we include an additional component 
to account for it in our spectral fitting (§ ||). 

For each galaxy, we extracted spectra in a number of 
concentric annuli, centred on the nominal X-ray centroid. 
We determined the centroid iteratively by placing a 0.5' 
radius aperture at the nominal galaxy position (obtained 
from NED) and computing the X-ray centroid within it. 
The aperture was moved to the newly-computed cen- 
troid, and the procedure repeated until the computed 



http: / /cxc. harvard.edu/ciao/thrcads/indcx. html 
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TABLE 1 
The galaxy sample 



Galaxy 


Type 


L B 


Lk 


Dist 


Scale 


Reff 


ObsID 


Date 


Exposure 






(1O 1O L ) 


(1O U L ) 


(Mpc) 


(" kpc- 1 ) 


(kpc) 




(dd/mm/yy) 


(ks) 


NGC 720 


E5 


3.1 


1.7 


25.7 


8.1 


3.1 


492 


12/10/00 


17 


NGC 1407 


E0 


6.4 


3.1 


26.8 


7.8 


4.4 


791 


16/08/00 


38 


NGC 4125 


E6 pec Liner 


4.7 


1.8 


22.2 


9.4 


3.3 


2071 


09/09/01 


63 


NGC 4261 


E2-3 Liner Sy3 


4.4 


2.2 


29.3 


7.1 


3.4 


834 


06/05/00 


34 


NGC 4472 


E2/S0(2) Sy2 


7.5 


3.2 


15.1 


11 


4.0 


321 


12/06/00 


34 


NGC 4649 


E2 


5.1 


2.5 


15.6 


13 


3.2 


785 


20/04/00 


21 


NGC 6482 


E Liner 


10.9 


3.2 


58.8 


3.6 


3.4 


3218 


20/05/02 


18 



The g^l^rin^.in the sample. Distances were obtained from ITonry 

do 



Note. 



zero-point (Jensen ct al 
Lb was obtained trom L 



2003) 
UK 



, except for NGC 6482, for which we adopted 
corrected to our distance. K B -b. md lumit i psiti 
= 5.48 and M K © = 3.41 (e.g. pVIaraston| [199f 



from 2MASS. We assumed M B 
the number of arc seconds corresponding to 1 kpc. We list the observation ID 
eliminated flaring intervals. 



ct al. (2001), corrected for the the new Ccphcid 
the kincmatical distance modulus from LEDA. 
es (Lk) and effective radii (R e ff) were obtained 
). We also list the image scale (Scale), which is 
(ObsID) and total exposure times, after having 



position converged. Typically the X-ray centroid agreed 
with that from NED. The widths of the annuli were 
chosen so as to contain approximately the same num- 
ber of background-subtracted photons and ensure there 
were sufficient photons in each to perform useful spectral- 
fitting. The data in the vicinity of any detected point- 
sources were excluded, as were the data from the vicin- 
ity of chip gaps, where the instrumental response may 
be uncertain. We extracted products from all active 
chips, excluding the S4, since it suffers from considerable 
"streaking" noise. Appropriate count-weighted spectral 
response matrices were generated for each annulus using 
the standard CIA tasks mkwarf and mkacisrmf. 

3.1. Background estimation 

One of the chief difficulties in performing spectral- 
fitting of diffuse emission is the proper treatment of the 
background. A set of standard blank-field "template" 
files are available for Chandra as part of the Caldb. We 
found, however, that the background template files are 
not sufficiently accurate to use in the very low surface 
brightness regions at large radii, which are crucial to de- 
termine interesting global mass constraints. The back- 
ground comprises cosmic, instrumental and non X-ray 
(particle) components. The cosmic component is known 
to vary from field to field, while the non X-ray back- 
ground exhibits long-term secular variability. To miti- 
gate the latter effect, several authors have adopted the 
practice of renormalizing the background template to en- 
sure good agreement with their data at high energies 
(> 10 keV). Such a procedure, however, also renormal- 
izes the (uncorrelated) cosmic X-ray background and in- 
strumental line features, which can lead to serious over- 
or under-subtraction. Given these reservations we chose 
to use an alternative background estimation procedure. 

Our method involved model ling the back g round , some- 



Buote ct al.| (|00|). All of 



what akin to the approach of 

the targets were centred on the ACIS-S3 chip, which is 
back-illuminated (BI) . To obtain constraints on the back- 
ground, we extracted spectra from a ~2' region centred 
on the SI chip, which is also BI, and from an annulus 
centred at the galaxy centroid and with an inner and 
outer radii typically ~2.5' and 3.3'. We excluded data 
from the vicinity of any point-sources found by the source 
detection algorithm. Although the diffuse emission from 



each galaxy typically had a very low surface-brightness 
on the SI CCD, we found that using two regions in this 
way with different contributions of source emission en- 
abled the background components to be most cleanly 
disentangled from the source. The A CIS focal plane also 
consists of front-illuminated (FI) chips, which have sig- 
nificantly different (and lower) background. To obtain 
an estimate of the background for these chips, we ex- 
tracted spectra from the entirety of each chip, excluding 
detected sources and data towards the edge of the chips 
where the exposure-map may be uncertain. 

In order to constrain the model, we fitted all spectra 
simultaneously, without background subtraction, using 
Xspec. Our model consisted of a single APEC plasma 
(to take account of the diffuse emission from the galaxy; 
the "source"), plus background components. These com- 
prised a power law with T — 1.41 (to account for the 
hard X-ray background), two APEC models with so- 
lar abundances and kT= 0.2 and 0.07 keV (to account 
for the soft X-ray background) and, to model the in- 
strumental and particle contributions, a broken power 
law model and two Gaussian lines with energies 1.7 and 
2.1 keV and negligible intrinsic widths. We have found 
that this model can be used to parameterize adequately 
the template background spectra. In general, the instru- 
mental contributions of the FI chips were very similar 
in shape. Therefore, the background components of all 
the FI chips were tied, assuming the normalization scaled 
with the spectral extraction area. For the BI chips, there 
was some evidence that the SI chip background can be 
somewhat larger at energies > 5 keV (although this is 
variable). In order to disentangle the source and back- 
ground components, given the general lack of photons 
in these spectra, we tied the abundances and temper- 
atures of the "source" APEC components between the 
extraction regions, but allowed the normalizations to be 
free. Where there was a significant improvement in the 
fit-statistic if this assumption was relaxed, we allowed 
the abundances or temperatures to fit freely. Notwith- 
standing, this assumption should not significantly affect 
our results. This model was able to fit all of the data 
well. In our subsequent spectral analysis, we did not 
background-subtract the data using the standard tem- 
plates, but took into account the background by using 
appropriately scaled versions of the models fitted to each 
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CCD, which were added according to the overlap between 
the source region and the CCD. We found that the stan- 
dard background templates fared much worse than these 
modelled background estimates when the data were from 
regions of low surface-brightness. We discuss the im pact 
of the background treatment on our results in 



We assess the impact of all these features in § 7.4 



7.3 



4. X-RAY IMAGES 

The X-ray image of each galaxy was examined to iden- 
tify any obvious surface-brightness disturbances or asym- 
metries which would be indicative of clear deviations 
from hydrostatic equilibrium. We note that low-level X- 
ray asy mmetries, such as th e "fingers of emission" identi- 
fied by Randall et al. (2004) in the adaptively-smoothcd 



images of JNCC 4649, probably do not merit concern 6 , as, 
provided care is taken to avoid seriously disturbed emis- 
sion regions, reliable mass p rofiles can be inferr ed even 
in mildly disturbed systems (Buote & Tsai 1995| ). 

In Fig we show the 0.1-10.0 keV 4 G7A-S3~Tmages of 
each of the systems. These images were first processed 
to remove point-sources, using the CI AO tool dmfilth, 
which replaces photons in the vicinity of each point- 
source with a locally-estimated background. NGC4261 
contains an AGN which appears as a bright central X- 
ray source and there is evidence of a small, low surface- 
brightness jet ( Zczas et al. 2005 ). We have also removed 
these sources from the image. The images were flat- 
fielded with the 1.7 keV monochromatic exposure-map 
(although this analysis is insensitive to the choice of en- 
ergy), and then smoothed by convolution with a 5" gaus- 
sian, to make large-scale structure more apparent. Due 
to the low surface-brightness nature of the emission at 
large radii, it is difficult to appreciate X-ray emission 
outside ~a few arc minutes in many of the images. How- 
ever, detailed spectral analysis and azimuthally-averaged 
surface brightness analysis reveals substantial hot gas ex- 
tending beyond the edge of the S3 chip in each system. 

None of the objects show very obvious disturbances 
in their X-ray emission o n the ACIS-S3 ch ip (such as 
those found in NGC4636: Jones et al. 2002). Some low- 
amplitude features are evident such as the faint jet in 
NGC4261 (which is not visible in the above images), a 
possible north-south asymmetry in NGC 1407 and some 
asymmetry, in particular an off- axis X-ray enhancement, 
in N GC 4125. Based on ad aptively-smoothed XMM im- 
ages, Groston et al. (2005) argued that the X-ray emis- 
sion in NGC 4261 is anti-correlated with the galaxy radio 
lobes. By inspection of the XMM images, this actually 
appears to be a very low-amplitude effect. It is not obvi- 
ous in the Chandra images, although the X-ray isophotes 
do align somewhat perpendicularly to the jet. In any 
case, this does not appear to have significantly disturbed 
hydrostatic equilibrium, since there is excellent agree- 
ment between our inferred mass profile and a model com- 
prising stellar plus DM components (§ |6|), which would be 
an extraordinary coincidence if hydrostatic equilibrium 
had been strongly disturbed. The limited field-of-view 
makes it difficult to assess asymmetries and disturbances 
on the other chips. NGC 4472 is known, however, to ex - 
hibit a disturbance outside ~6' fllrwin fe Sarazin 1996 ). 



5. SPECTRAL ANALYSIS 

Spectral-fitting was carried out in the energy-band 0.5- 
7.0 keV, to avoid calibration uncertainties at lower ener- 
gies and to minimize the instrumental background, which 
dominates at high energies. The spectra were rebinned 
to ensure a S/N ratio of at least 3 and a minimum of 
20 photons per bin (to validate x 2 fitting). We fitted 
data from all annuli simultaneously using Xspec. To 
model the hot gas we adopted a vapec component, plus 
a bremsstrahlung component for all annuli within the 
twenty-fifth magnitude isophote (D25) of each galaxy, 
taken fro m the Third Reference Cata log of Bright Galax- 
ies (RC3: ie Vaucouleurs et al. 1991), to account for un- 
detected point-sources (this model gives a good fit to the 
composit e spectrum of th e detected sources in nearby 
galaxies: [rwin et al. |2003 ). We used a slightly modified 
form of the existing Xspec vapec implementation so that 
Zp e is determined directly, but for the remaining elements 
the abundance ratios (in solar units) were directly deter- 
mined with respect to Fe. This was useful since, in gen- 
eral, the data did not enable us to determine any abun- 
dance ratio gradients and so we tied the abundance ratios 
between all annuli. Where abundances or abundance ra- 
tios could not be constrained, they were fixed at the Solar 
value. We adopted the sol a r pho tosphcric abundances 



standard of Asplund et al. ( 2004 ). We r efer the inter 
ested reader to [Humphrey fe Buotc ( 2006 ) for a detailed 
discussion of this choice and how to convert our results to 
older abundance standards. In the interests of physically 
reasonable results, we constrained all abundances and 
abundance ratios to the range 0.0-5.0 times solar. The 
absor bing column density (Ah) was fixed at the Galactic 
value ( Dickey fe Lockman 199C ) ; the effect of varying Ah 
is discussed in § |7.6| . For NGC 4261, our innermost an- 
nulus contained substantial contamination from the cen- 
tral AGN. However, this was sufficiently absorbed that 
the thermal emission from the gas can be clearly disen- 
tangled from it. To account for the AGN emission, we 
fitted a highly absorbed (ATh= 10^4 X 10 22 cm~ 2 ) power 
law component (r — 1.4 ± 0.8). We dis cuss the impact 
of including this annulus on our fits in § \IA. 



but mass analysis inside this radius should be reliable. 

6 Although the authors suggested these may arise from bulk 
convective flow, the spectra do not agree with simulations of such. 



To account for projection effects, we used the pro jet 
model implemented in Xspec. This model, unfortunately, 
does not take into account the emission from gas out- 
side the outermost shell, which is also projected into the 
line-of-sight. To take account of this effect, we assumed 
that the emission outside this shell has the same spectral 
shape as the emission in that shell and a dens ity profile 
well-described by a /3-model (e.g. Buotc] 2000| ). We in- 
cluded an extra spectral component to our fits of each 
annulus to account for projected emission from this gas. 
To estimate the parameters of the /3-profile, we we fitted 
the galaxy surface brightness, using dedicated software, 
in the 0.1-3.0 keV band. Although a single /?-model did 
not always match the fine detail of the surface bright- 
ness profiles, it adequately parameterized the data for 
our purposes (our results are not expected to be strongly 
dependent upon the parameters of this fit). 

We obtained good fits to the spectra of each galaxy 
with this model. The best-fitting abundances were in 
exc ellent agreement with thos e of other early-type galax- 
ies (Humphrey & Buotc [2006 ), and are shown in Table ||. 
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Fig. 1. — X-ray images of each of the galaxies in the sample. None of the systems show evidence of large-amplitude disturbances which 
would indicate a. violation of hydrostatic equilibrium. Some lower-amplitude asymmetries do persist in some of the images, which we discuss 
in detail in § U^4. 



We note that |Randall ct al.| ( |2006D found Z sl /Z Fe ~1.7 for 
NGC 4649 (adjusting to our abundances standard) when 
fitting the data from single, large aperture, which they 
argued points to substantial enrichment from SN II, in 
stark contrast to the pr edominantly SN la enrichm ent we 
found in such galaxies ( Humphrey fc Buotc 200€| ) . From 
our analysis, however, Zsi/Zp c ~l, which is more consis- 
tent with our results for other systems. The discrep- 
ancy appears to be related to the "Fe bias" (where Zp e 
is systematically underestimate d if one assu mes multi- 
temperature gas is isothermal: Buotc 2000[ ) which has 
suppressed their large aperture Zp c in comparison to 
their spatially-resolved results (which agree better with 
our measurement). 

Error-bars were computed via the Monte-Carlo tech- 
niqu e which we have e xtensively used in past analyses 



(e.g. Buote et al. 2005 ). We simulated spectra from the 
best-fit models, which were then fitted exactly analo- 
gously to the real data. We performed 25 simulations, 



which were sufficient to assess the distribution of the fit 
parameters about the best-fit values; the standard de- 
viation of this distribution corresponds to the 1-er con- 
fidence region. Assuming that we have found the global 
minimum, and the fit statistic follows a y 2 distribution 
this is statistically equivalent to searching the parameter 
space for changes in the fit statistic. Tem perat ure and 
density profiles are discussed below (§ 6.2 and § 6.3) 

6. MASS MODELLING 

6.1. Assumed potential method 

We adopted two complementary approaches in order to 
determine the mass profiles of the galaxies in the sam- 
ple. The first method, discussed here, was found to be 
less sensitive to the assumptions of the modelling and 
therefore was adopted as ou r default. We discuss our 
alternative approach in § p. 4. 



Starting with a parameterised model for the tempera- 
ture (T) and gravitating mass (M grav ) profiles, the equa- 
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TABLE 2 

Emission-weighted average abundances 



Galaxy 



X 2 /dof 



Zfo 



Zq/Zi 



Zno/Zfc ZMg/ZF c Zsi/Zi 



Zs/Zj 



Zni/Zf 



NGC 720 1 


383.4/357 


0.80±° ;« 


0.30+.0.28 


0.68±0.67 


1.26±0.35 








NGC 1407 1 


222/221 




37+ - 21 

u " 3 '-0.25 




1.10±0.23 


i 91 +0.31 
-0.27 


2.2±1.1 


o 3+1.7 
°- -1.3 


NGC 4125 


327/307 


55+ - 22 
u - oo -0.13 


29+ 13 
u - za -0.09 


0.62±0.14 


0.33±0.20 








NGC 4261 


307/319 


1.72±0.50t 


<0.23 


36+ 79 


0.83±0.23 


1.2±0.4 




1 8+ 2 ' 3 


NGC 4649 


563/491 


r, , 9 + 0.87 

^••"-0.37 


< 0.15 




0.97±0.13 


1.02±0.13 




1 42+ ' 85 


NGC 4472 1 


785/740 


1 -*-0.4 


0.51+.0.12 


0.95 ±0.44 


1.02±0.11 


1.25+.0.11 


2.36±0.33 


3.28+.0.61 


NGC 6482 


256/262 


>2.5 


0.34+.0.20 




1.15±0.18 


1.3+.0.3 




n + 1-5 

2 



Note. — The best-fitting globally-averaged emission-weighted abundances and abundance ratios for each galaxy, shown 
along with the quality of fit. Statistical errors represent the 90% confidence region. Where we were able to constrain an, 



ahnnHanrn gradient, we estimated an emission-weighted Zfc. cxtrap 
2006); those affected galaxies arc marked 1 — results taken from 

not be constrained, they were fixed at the Solar value, and listed as 



erture (sec Humphrey 8z Buote 



Humphrey &: Buotc (2006). Where parameters could 
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Fig. 2. — Deprojected temperature profiles of each galaxy, ordered by M v i r (Table ^). The data-points are shown, along with the best-fit 
parameterized model determined from simultaneously fitting the p g and temperature profiles (see text). Where data-points were excluded 
from the fitting, they are denoted by dashed lines. Errors shown are 1-cr. 



tion of hydrostatic equilibrium can be solved for p g thus: 

In ( £*-) = - in (pi - G,m p [ R M "£< V d R 
\P 9 oJ \T J * p J Ro kTR 2 

where R is the radius from the centre of the gravita- 
tional potential, p g is the gas density, p a o and Tq are 
density and temperature at some "reference" radius Rq, 



k is Boltzmann's constant, G is the universal gravita- 
tional constant, m p is the atomic mass unit and /i is the 
mean atomic weight of the gas. In our fitting we explic- 
itly ignored the contribution of the gas to the gravitat- 
ing mass, but we subsequently verified this contributed 
< 1% of the total gravitating matter within 100 kpc, jus- 
tifying this assumption. We developed software to fit p g 
and temperature profiles simultaneously using this proce- 
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dure. For speed we assumed that the density and temper- 
ature data-points were each evaluated at a single point, 
the radius of which was given by: 

Ri = (0.5 * (Rinj- 5 + Rout] 5 )) 2/3 (2) 
where Ririi an d Routj were the inner and outer radius 
of the bin (see Lewis et al. 1 12003|) - 



6.2. Temperature profiles 

There were considerable differences in the temperature 
profiles from object to object (Fig|), so that we were not 
able to adopt a "universal" profile for all of the systems. 
A priori we do not expect any particular form for the 
temperature profile and so we determined appropriate 
functional forms for our temperature models empirically. 
Based on experience, the following "toolbox" of models 
provided adequate flexibility to ensure at least one model 
can describe the tem perature profiles reasonably well (sec 



Buote et al. 2006b) 



T = T + T 1 [1 + 
T=[T +T lX Pl }< 



T-- 



A 



A + B 

B 
A + B 



T 2 x P2 1-e 



l + Xi 

[T 2 + T 3 (l + x 2 )- p *] 



(3) 
(4) 



(5) 



where x = (r/r c ), x x = (r/r cl ), x 2 = (r/r c2 ), A = 
(1 + r/r tl )- 3 ^ and B = e(l + r/r t2 )-^\ T , T x , T 2 , 
?3, r c , r&, r c2 , m, r t2 , pi, p 2 , Pe and e are parameters 
of the fit. For NGC4261, we ignored temperature data- 
points from 15-25 kpc, which were poorly-determined 
and seemed erroneous. We experimented with fitting the 
projected (rather than deprojected) spectra, and found 
no evidence of any features (in either temperature or den- 
sity) around this range of radii, strongly implying that 
they arise solely due to deprojection noise. The temper- 
ature profiles and best-fit models are shown in Fig g. 

Our deprojected temperature profiles generally agree 
with those appearing in the literature for these objects. 
(Although most of these are projected profiles, typi- 
cally deprojection does not strongly alter the overall 
charac ter of the temperature profile.) p' Sullivan et al 



(2003) reported ROSAT profiles for all of the galaxies 
which, although substantially less well-constrained, agree 
well with our results. Likewise our NGC4472 tempera- 
ture profile agrees well with the (les s well- constrained) 
ROSAT profile of |frwin fc Saraainl ( |f~996| ). Our pro- 



file for NGC 4649 is in reasonabl e agreement w i th th e 



projected XMM measurements of Randall et al. ( 2006] ), 
and likewise our measured p rofile of NGC 6482 agree s 
with the deprojected results of Khosroshahi et al. ( 2004| ). 
Our temperature profiles for NGC f 407, NGC 720 and 
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N GC 4472 were also in a greem ent with those we reported 
in 



Humphrey & Buotc ([2006 ) 



6.3. Mass-fitting results 

We tested three different mass-models against the data. 
In order to investigate the suggestion that historically 
large c values found based on X-ray analysis were an 
artefact of the omission of t he stellar mass, as we ll as 
to investigate the scenario of Loeb & Peebles (2003), we 
first tested a model comprising a single NFW profile. Al- 
tho ugh stellar kinem atical results would seem to rule out 



the Loeb fc Peebles picture, our analysis of more mas 



sive systems (Gastaldello et al. 2006) does suggest that 
the stellar mass may not be uniformly required in all 
systems. In order to take into account the stellar mass, 
we fitt ed a mode l comp rising an NFW DM component, 
plus a Hcrnquisl] ( |199C , hereafter H90) mass component, 
the R c ff of which being fixed to that measured in the 
K s -band (Table @). The H90 model is, in projection, a 
good approximation to the familiar de Vaucouleurs pro- 
file of elliptical galaxies. To test whether the DM halos 
retained any evidence of their response to baryonic con- 
densation, we further adopted an H90 component, plus 
an NFW com ponent modifi e d by the adiabatic contrac- 



tion model of Gnedin et al. (2004) 7 . Hereafter, we refer 



to these three models as, respectively, NFW, NFW+stars 
and AC NFW+stars. Our computed M v j r for each sys- 
tem included both dark and stellar mass. For NGC 4472 
and NGC 4649, which lie in Virgo, there is the possibil- 
ity that the DM halo may have experienced some tidal 
truncation at a radius <R v i r - Our measured Virial quan- 
tities relate to the original halo prior to truncation. We 
also experimented with replacing t he N FW component 
with the less cuspy Navarro et al. ( 2004 , hereafter N04) 
model, which gives an improved fit to DM halos in high- 
resolution N-body simulations. However, since the M v j r -c 
relation was calibrated using the NFW model we treat 
t his c hoice as a systematic effect and it is discussed in 



7.1. For NGC 4261 we ignored a deviant p g data-point 
at ~11 kpc, in addition to the excluded temperature 
data-points discussed above. 

In Fig H we show the density profiles (along with the 
best-fitting AC NFW+stars model, which is described 
below). In Fig|^(a) we show the best-fitting 1-er contours 
of c versus M v i r for the NFW model fitted to each galaxy. 
The fits were typically, but not uniformly, poor (Table ||). 
We found very large (;§>20) values for c, completely in- 
consistent with the expectation of N-body simulations. 

Good constraints on the global halo properties typ- 
ically require interesting density and temperature con- 
straints over as large a radial range as possible. In our 
case, the absence of data outside ~50-100 kpc (~0.1- 
0.2R V j r ) therefore makes the inner data-points critical 
in determining the profile of the halo. Unfortunately, 
since the scale radius of a galaxy-size DM halo is ~10- 
30 kpc, there is some degeneracy between the DM and 
s tella r mass components at small radii. As we discuss in 



§6.5 there are considerable uncertainties in estimating a 
reliable mass-to-light (M/L) ratio from the characteris- 
tics of the stellar population. We found that the results 
are extremely sensitive to the stellar M/L adopted; we 



found that varying this ratio by as lit tle as 20% could 
cause M v ; r variations of ~50-100% (see Humphrey et al 
2005 ). It was therefore necessary to allow the stellai 



mass to be determined as a parameter of the fit. This, 
unfortunately, made it very difficult to constrain M v i r or 
c, unless additional constraints were applied to the fit. 

One way to achieve this is to constrain the fit to lie 
on the mean M v j r -c relation predi cted from N-body sim- 



ulations (e.g. Bullock et al. 2001). Although this would 
prevent our measuring M v ; r and c independently, it would 
enable us to determine whether the data are consistent 
with the mean relation. However, this relation was de- 
termined for an unbiased sample of DM halos, whereas 
our selection criteria (§ ||) should bias us towards sys- 
tems which have not recently had a merger (implying 
earlier- forming, hence more concentrated, objects). Fur- 
thermore, individual halos are not expected to lie exactly 
on the mean M v j r -c relation, but be scattered about it. 
Nevertheless, we experimented with applying this con- 
straint. The data for each galaxy were consistent with 
this model, but we found M v i r was generally poorly con- 
strained, and extremely sensitive to any scatter we intro- 
duced about the mean M v ; r -c relation. 

A far more useful way to constrain the fit was to re- 
strict the total baryon fraction (fb) in the system. Such 
a constraint is useful since we found that, for a given 
system, fb determined from our fits was strongly anti- 
correlated with the measured M v ; r . To estimate fb, we 
computed the gas mass by extrapolating our p g model 
from the centre of the innermost radial bin to the Virial 
radius. The contribution of stars to the total baryon frac- 
tion was derived from the stellar mass found by our fits. 
We crudely took into account the fact that not all of the 
stellar mass within R v ; r is necessarily contained in the 
central galaxy by scaling this mass by the ratio of the to- 
tal B-band light of all pu t ative "group" members listed in 
the catalogue of Garcia ( 1993| , hereafter G93) to that of 
the central galaxy. This is likely to overestimate slightly 
the stellar mass content, since it assumes the same stellar 
M/L ratio for all low-mass companions/ group members, 
whereas some fraction of these are likely to have substan- 
tial young stellar populations, with lower M/ L r atios. 
We discuss the impact of this assumption in § G93 
lists NGC 4649 as belonging to the NGC 4472 Voup" , 
whereas they both have distinct X-ray halos, indicating 
they are, in fact, distinct systems. As a zeroth order ap- 
proximation, we therefore assumed that the total B-band 
luminosity was divided between the two "subgroups" in 
proportion to the central galaxy's B-band luminosity. In 
practice, between ~25% (for NGC 4472) and 84% (for 
NGC 4125) of the B-band light of the system resides in 
the central galaxy. NGC 6482 was not listed in G93, but 



as it is known to be relatively isolated ( Khosroshahi et al 
2004), we assumed that ~80% of its mass is in the cen- 



7 Available publicly from http://www.astronomy.ohio- 
state.edu/~ognedin/contra/ 



tral galaxy, consistent with the other relatively isolated 
systems. 

Based on hydrodynamical simulati ons incorp o rating 
gas cooling and supernovae feedback, Kay et al. (2003) 
predicted fb as a function of Virial temperature for sys- 
tems with M v i r >a few xl0 12 M Q . Fitting their data by 
eye, converting from Virial temperature to mass, and as- 
suming a Universal baryon fraction of 0.16, we estimate 



ft = 0.062 log 10 (M 14 ) + 0.13 (M u < 1.02) (6) 
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Fig. 4. — 1-cr (M v i r ,c) confidence regions for each object. Panel (a) shows the confidence contours found when fitting the NFW model. 
In panel (b) we show the contours for the AC NFW+stars model, restricting fb by Eq (see text), and panel (c) shows contours for the 
same model, but now restricting 0.032<f] 3 <0.16. The results for the NFWj±staES mode] wpr> .very similar to those for the AC NFW+stars 
model. On each plot, we also show the prediction from the "toy model" of 
in the simulated DM halos (the region bounded by the dotted lines). 



Bullock ct alJ (2001) (solid line) and the approximate 1-cr scatter 



TABLE 3 
Quality of the mass fits 



Galaxy 


NFW 


NFW+stars 


AC NFW+stars 


NGC 720 


1.9/9 


1.0/8 


0.9/8 


NGC 1407 


26.7/9 


20.7/8 


20.7/8 


NGC 4125 


23.4/11 


9.5/10 


10.8/10 


NGC 4261 


22.6/12 


14.0/11 


14.0/11 


NGC 4472 


35.2/20 


34.9/20 


35.2/20 


NGC 4649 


30.2/7 


11.0/6 


11.5/6 


NGC 6482 


0.5/5 


2.4/4 


1.7/4 


Note. — 


The x 2 /dof of the fits 


to the density and 



temperature profiles used to infer the mass, for the three 
basic mass-models adopted. For the NFW.4- stars and AC 
NFW+stars models, we constrain f^ to Eq la. 



0.14 



(M u > 1-02) 



where M14 = M v i r /10 M@, with an approximate scat- 
ter of ±0.02. Adopting this constraint (including the al- 
lowed range of scatter) and fitting the NFW+stars model 
resulted in significant improvements in the fit quality 
over the simple NFW model (Table ||). We show the 
resulting c-versus-M vlI contours in Fig |(b), and sum- 
marise our results in Table [|. Clearly adding the stellar 
mass component allows the DM halos to be substantially 
less concentrated, since less DM is required in the centre 
of the halo. These results were in much better agreement 



with the results of N-body simulations than those ob- 
tained with the NFW model. There is, however, a slight 
trend towards more concentrated halos than ACDM. 

Since the results of simulations can be sensitive to 
the rather uncertain process of feedback, we addition- 
ally adopted, as a somewhat less restrictive constraint 
on fb, 0.03<fh< 0.16, the l ower limit being ^the lowest 
values found in Kay et al.'s simulations. The shallower 
potential of very low-mass halos makes it more difficult 
for them to hold onto their hot gas, and so our lower 
limit on fb may be an overestimate if the Virial mass is 
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TABLE 4 
Best-fitting NFW+stars results 



Galaxy 


f b =Eq H 








0.032<f b <0.16 












M vir (lO 42 M ) 


R vir (kpc) 


c 


fb 


M vir (10 12 


M @ ) 


R vir (kpc) 


c 


fb 


NGC720 


6 - 6 -3.0 


480l|° 


18.+ 30 - 


n n44+ 037 
u.u44_ 003 


6.C 


+ 6.0 
-4.3 


480 ± 120 


18. 


+49. 
-10. 


n 044+ ' 095 
u.u44_ 012 


NGC1407 


16. ±6. 


fjrcn+80 
DOU -100 


18+i L 


065+ - 041 


21. ± 15. 


720+ 140 
,zu -200 


15. 


+ 16. 
-6. 


032+ 130 


NGC4125 


6.2±°;f 


470+ 20 , 


10. j^; 


oq+0.035 


7 2 +1.4 


500+ 30 
OUU -160 


9.3 


+ 11. 

-2.1 


032 +013 


NGC4261 




1040±|°° 


3.7 ±1.7 


14+ 001 


57. 


1 260 
-15. 


oqo+760 
9MU -100 


4.0 ±2.0 


1fi+ 000 
u - iD -0.13 


NGC4472 


QO +6. 

oo --10. 


«*>±iao 


13.+ 4 ' 


0.084±g:8H 


63 


+ 17. 
-44. 


1020+ 90 
iuzu_ 300 


10. 




0.032t°J 3 


NGC4649 


?1 +7- 
°°--13. 


840if° 20 


21+| : 


0.086±g;ggl 


93 


+26. 
-73. 




15. 


+11. 

-3. 


0.032«J2 


NGC6482 


7 1+ 4 ' 4 
— 1. V 


50011° 


18+J?" 


0.075l°;° 3 3 2 


3.C 


+5.5 
-1.5 


390^° 


38. 


+76. 
-24. 


lfi+ - 00 



Note. — The best-fitting results for the NFW+stars model. All error-bars shown correspond to 90% confidence regions. The 
fit results for the AC NFW+stars model are very similar, and arc shown in Fig hi Results are shown for the fits using the two 
different constraints on f^ we adopted (sec text). 
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Fig. 5. — The K-band mass-to-light profile of each galaxy. Results are shown both for the NFW+stars and the AC NFW±stars models. 



For the latter the shaded region indicat 
profile" mass-modelling described in 
independently. 



^s the approximate 1-cr confidence region. We also show data-points derived from the "parameterized 



6.4. We stress that these data-points are not fitted by the models shown here, but they are derived 



small. However, imposing a lower limit on fj, in the fitting 
algorithm actually works to exclude the most massive so- 
lutions, which we would exp ect to be closer to baryonic 
closure ( Mathews et al. 2005 ) . It is conceivable that some 
more massive (M v i r > 10 iJ Mo) systems are rare exam- 
ples of "dark" groups which have unusua lly l ow baryon 
fractions, a possibility we return to in S_|S.7| . Notwith- 
standing, the results are shown in Figs jj(c), which are 



qualitatively similar to those obtained applying the more 
restrictive constraint on fb, although with larger uncer- 
tainty. 

In Fig H we show the gravitating mass to K-band light 
(M grav /LK) ratio profile implied by our best-fit models 
for each system. In each case we found that the X-ray 
emission was considerably more extended than the opti- 
cal light. We also show data-points estimated from "pa- 
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rameterized profile" mass modelling (§ 3.4), which tend 
to agree reasonably well; the slight systematic differences 
between the profiles are an artefact of the assumptions 
used to der ive the data-points and we discuss this in de- 
tail in § 3.4 Clearly M grav /LK increases very slowly with 
radius within R c ff , rising very steeply outside this range. 
This arises naturally from the very different shapes of 
the stellar and DM halos, and is similar to M/ L profiles 
seen from st el lar ki nematics and the results of jBrighcnti 
& M jathewsj fli~997| ) for NGC4472 and NGC4649. By 
Rvir, Mg rav /LK reaches as high as ~20-40 M Q /L Q for the 
galaxy-scale systems or ^100-200 M Q /L© for the group- 
like objects. We stress that this only includes the light of 
the central galaxy which, for the group-like systems may 
be a little as ~25% of the total luminosity. 

6.4. Parameterized profile mass modelling 

We briefly discuss here an alternative technique to de- 
termine the mass profiles of X-ray bright objects which 
we have extensively employed in our previous studies, as 
well as the companion papers to this pre sent work (e.g . 



Lewis ct al. 


2003; 


Buotc et al. 


2006b 


2006 




Zappacosta et al. 


2006) 


. This 



we here dub the "parameterized profile" method involves 
parameterizing independently the temperature and den- 
sity profiles of the system with simple, empirical models. 
These functions were then inserted into the equation of 
hydrostatic equilibrium, which we solved for the mass en- 
closed within any given radius. The temperature pro files 
were parameterized with the models discussed in § 3.2, 
and to fit p g we adopted, where appropriate a, /3-model, 
a "double-/3" model or a "cusped-/3" model, defined, re- 
spectively, as: 



-3/3/2 



Pg = PgO [l + {r/r c ) 2 

P g = \/pg° [1 + (r/r c ) 2 r 30 + Pgi [1 + (r/r c2 ) 
Pg = p a0 2^-^(r/r c r* [l + ( r /r c f]- 3f3/2+e/2 



(7) 

-i : '' 02 (8) 
(9) 



Where the parameters p g o, p g ±, r c , r C 2, /3, @2 and e are 
determined by the fit. Fitting these models to the sim- 
ulated temperature and density profiles discussed in § || 
(which were used therein to estimate the error-bars on 
kT and p g in each data-bin) allowed us to estimate the 
scatter in the mass data-points arising from statistical 
noise, and hence the error-bars. 
For a full discussi on of this te c hnique , we refer the 



interested reader to Buote et al. (2006b), who demon 



strate the good agreement typically found between this 



method and the assumed potential modelling of § 3.1 
when fitting high-quality data. However, the mass data- 
points, especially at the innermost and outermost radii, 
are rather sensitive to the parameterized models adopted 
to fit the temperature and, especially, p g . The system- 
atic uncertainty introduced by the choice of p g model can 
be considerably larger than the statistical error. For our 
purposes the absence of data at very large radii, which 
are vital to constrain the curvature of the mass model, 
exacerbated by the uncertainty introduced at small radii 
due to the uncertain stellar mass-to- light ratio, magnified 
the impact of these systematic effects. Notwithstanding 
these reservations, it is still interesting to compare the re- 
sults obtained via both mass-fitting methods. We show 



m Fig | the mass data-points computed using parame- 
terized potential mode lling , along with the best-fitting 
mass models found in § 6+ . Clearly there is good overall 
agreement between the two methods although there are 
some systematic differences, which reflect the systemat- 
ics inherent in our choice of parameterized model for p g . 

6.5. Stellar mass-to- Light ratios 

It is interesting to compare the stellar M/L ratios 
(M*/L) determined by our fitting to the expectations of 
stellar population synthesis models. In order to ensure 
that the optical light traces the stellar mass as closely as 
possible, we opted to perform this comparison in the K- 
band. Table shows M*/Lk determined from our models 
using eq ^ to constrain ft,. Since AC tends to increase 
the cuspiness of the DM profile we found a significantly 
lower mass-to-light ratio for the AC NFW+stars model 
than for the NFW+stars model. 

To compare our measured M»/Lk to single burst stel- 
lar population synthesis predictions, we first estimated 
a mean emission-weighted stellar age and metallicity for 
each galaxy, as outlined in Appendix |A|. We linearly 
interpolated synthetic M ^/Lk valu e s bas ed on the stel- 



lar population models of |Maraston ( 1998[ ) from updated 
model-grids made available by the author 8 . For typ- 
ical early-type galaxies, K-band and 2MASS K s -band 



magnitud es should differ by <0.1 magnitudes ( [Carpen- 



ter 



2001), so we were able to compare directly the syn 
thetic K-band M/L ratios with our measured M*/Lk ra- 
tios. The predicted M*/Lk ratios are shown in Table || 
for different assumptions about the stellar IMF, which is 
poorly-known in early- type galaxies. In this case we show 
predicted M*/Lk a ssuming a stan dard Salpeter IMF, 
and for the IMF of Kroupa (2001). It is immediately 
clear that these ratios are very sensitive to this choice; 
M*/L K is typically ~50-60% higher if the Salpeter IMF 
is used. 

Our measured M*/Lk for the NFW+stars models are 
typically ^20% lower than the synthetic M/L ratios, 
assuming the Kroupa IMF. Using the AC NFW+stars 
models, the discrepancy is ~40%. Assuming a Salpeter 
IMF, the discrepancies for both models are consider- 
ably larger. This would seem to rule o ut the Salpeter 
IMF, i n agr eement with the conclusions of Padmanabhan 



t al, ( 2004 ). The best-fitting M v ; r and c are sensitive to 
M*/Lk- If we fix M*/Lk to the synthetic value, this es- 
sentially pushes all the galaxies, except NGC4261 (for 
which the measured and synthetic values are in excellent 
agreement) and NGC 720 in the direction of the high- 
M v ; r range of their confidence contours shown in Fig ^. 
For NGC 720, M v j r is lowered and c increased. The fits 
are then typically much worse (A% 2 ~7-35), and the loci 
in the M v ; r -c plane slightly more discrepant with simula- 
tions. 

There are a considerable number of systematic uncer- 
tainties in the computation of the synthetic M/L ra- 
tios, not the least of which is the very uncertain IMF, 
which could probably account for the mo des t discrep- 
ancy with our NFW+stars results (see § 8.4). In the 
case of NGC 720, the rather young age inferred for the 
stellar population (~3 Gyr) leads to a significantly lower 



' Ihttp:/ /www-astro. physics. ox. ac.uk/~maraston/Claudia' 
effiT - 



_SteIlar_Population_Modcls.html 
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Fig. 6. — Mass profiles for each galaxy. The data-points were computed using parameterized profile modelling (§ |6.4j ). In addition we 
show the best-fit NFW+stars mass models from assumed potential modelling, which generally agree reasonably well with the data-points, 
indicating the consistency of both approaches to determine the mass profiles. We show the total gravitating mass model (solid line) and, 
separately, the stellar mass contribution (dotted line), the DM contribution (dashed line) and the gas mass (dash-dot-dot line). The models 
are extrapolated out to R v i r . Errors shown are l-cr. 



TABLE 5 
Stellar mass-to-light ratios 



Galaxy Lk/Lb 



Fitted M*/L K (Mq/Lq) Pop. synthesis M*/L K (M /L ) 
NFW+stars AC NFW+stars Salpeter IMF Kroupa IMF 



NGC 720 
NGC 1407 
NGC 4125 
NGC 4261 
NGC 4472 
NGC 4649 
NGC 6482 



5.5 
4.8 
3.8 
5.0 
4.3 
1.9 
2.9 



7+0.52 
-0.71 
H0.25 
-0.32 

0.72+0.11 
1.2+0.1 
51 + 24 

0.90+0.1.3 
F0.18 
-0.27 



0.77] 

0.52^ 



0.73] 



0.54- 1 



,+0.42 

-0.48 

0.35+0.25 

0.53+0.11 

1.0+0.1 

0.36+0.1 

0.65+0.12 

,+0.19 



0.52 i 



0.23 



0.54+0.11 
1.6+0.2 
1.7+0.5 
1.9+0.1 
1.3+0.3 
1.7+0.2 

1.58+0.02 



0.35+0.07 
1.1+0.1 
1.1+0.4 
1.3+0.1 

0.83+0.15 
1.1+0.1 

1.05+0.02 



Note. — K-band stellar mass-to-light ratios measured from our fits to the data using both 
the NFW+stars and the AC NFW+stars models. Since AC tends to increase the cuspiness of 
the DM profiles, M*/Lk is substantially lower for the AC NFW+stars models. We also show the 
predicted Mi/ t -i; y a , 1l1pg .derived from simple stellar population synthesis, assuming cither the 
Salpeter or |Kroupa (2001) IMF. 



synthetic M»/Lk than measured. Fitting temp late mod 



els to s patia lly- resolved spectra of this system, Rembold 
et al| (|2005| ) found evidence of a significant age gradi- 
ent, which falls from ~12 Gyr in the centre to ~3 Gyr 
by 1 kpc. This may, therefore, represent a system in 
which a relatively small fraction of the stellar compo- 



nent, produced in a modest, recent star- formation event 
("frosting") dominates the optical line emission. In this 
case, the synthetic M,/Lk may be underestimated. We 
return to this issue in § p.4| . 

7. SYSTEMATIC ERRORS 
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TABLE 6 
Systematic error budget 



Galaxy 




Astat 




Astars 


A KV*1 


Aasym 




Aspcctra 


Adist 


Abaryons 












M vir /10 12 


M 










NGC720 


6.6 


+2.4 
-3.0 


+0.06 


+0.5 

-0.3 


-1.7 


-0.2 


-0.9 


+U.7 
-0.2 


+0.9 
-0.6 


-0.01 


NGC1407 


16 


±6 


-0.1 


+ 1 
-0.5 


-3 


+3 


+i 

-3 


+ 1 
-3 


+3 
-4 


-0.3 


NGC4125 


6.2 


+0.8 
-2.3 


-0.2 


+ 1-7 
-2.1 


+0.09 
-0.2 


+0.2 
-0.03 


-1.0 


+0.5 
-2.0 


+0.7 
-0.6 


-0.3 


NGC4261 


67 


+41 
-15 


+2 


+ 15 

-24 


+9 


+5 
-0.5 


+0.08 
-0.6 


+ 13 
-0.04 


+4 
-7 


-3 


NGC4472 


33 


+6 
-10 


+0.5 


+6 

-0.05 


+2 


+12 


-4 


+8 
-1 


±3 


-5 


NGC4649 


35 


+7 
-13 


-8 


+ 16 
-3 


+2 


-1 


+63 
-2 


+6 
-2 


±3 


-19 


NGC6482 


7.1 


+4.4 
-1.7 


-0.8 


+ 1.9 
-0.3 


-0.1 


+0.8 


+0.1 
-3.6 


+0.3 
-0.8 


+0.9 
-0.8 


-0.2 


c 


NGC720 


18 


+3U 
-8 


+0.1 


±2 


+6 


+0.2 


+2 


+2 
-1 


±2 


+0.03 


NGC1407 


18 


+ 11 

-7 


-0.8 


+2 
-3 


+4 
-3 


-4 


+4 
-1 


+6 
-3 


+5 
-3 


+0.02 


NGC4125 


10 


+5 
-2 


-1 


+5 
-4 


+0.9 
-0.1 


+0.05 
-1.0 


+2 


+2. 
-0.7 


+ 2 


+0.3 


NGC4261 


.3.7 


±1.7 


-1.4 


+3.6 
-1.2 


+ 0.10 
-0.7 


-0.3 


+0.08 
-0.03 


-1.5 


+0.8 
-0.1 


+0.08 


NGC4472 


13 


+4 
-2 


-2 


+0.4 
-2 


-0.7 


-5 


+ 1 


+0.2 
-3. 


±2. 


+1.0 


NGC4649 


21 


+<; 

-3 


-2 


+3 
-6 


-0.9 


+2 


+3 
-4 


+0.4 
-3 


±3 


+7 


NGC6482 


18 


+ 13 

-8 


+ 1 


+2 
-6 


+3 


+6 


+51. 
-0.4 


+2 


±3 


-0.7 



+0.52 
-0.71 
+0.25 
-0.32 
±0.11 

±0.1 

+0.24 
-0.31 

±0.13 

+0.18 
-0.27 



M. 
+0.28 

-0.18 
+0.71 
-0.16 
+0.62 
-0.20 
+0.7 
-0.8 
+0.45 
-0.06 
+0.80 
-0.19 
+ 0.50 

-0.04 



NGC720 

NGC1407 

NGC4125 

NGC4261 

NGC4472 

NGC4649 

NGC6482 



0.77 

0.52 

0.72 

1.2 

0.51 

0.91 

0.73 



+0.17 

+0.06 

+0.04 

+0.05 

+0.06 

+47 

+0.13 



-0.15 

+0.12 

-0.04 

+0.01 

-0.03 

+0.008 

-0.001 

+0.05 



-0.06 

M«/L K (AC NFW+stars 



+0.07 


-0.04 


+0.14 


-0.19 


+0.06 


-0.04 


+0.04 


+0.0009 


-0.05 


-0.006 


+0.20 


-0.01 


-0.06 


-0.16 


-0.15 


-0.51 



±0.10 

+0.09 
-0.18 

+0.05 

+0.09 

-0.010 

+0.13 

-0.02 

+0.05 

-0.009 

-0.05 



±0.16 

+0.10 
-0.08 

±0.15 
±0.2 

+0.11 
-0.09 

±0.18 
±0.15 



-0.001 

+0.001 

-0.003 

-0.006 

-0.01 

-0.03 

-0.01 



NGC720 


0.54 


+0.42 
-0.48 


-0.12 


NGC1407 


0.35 


±0.25 


-0.06 


NGC4125 


0.53 


±0.11 


-0.04 


NGC4261 


1.0 


±0.1 


+0.02 


NGC4472 


0.36 


±0.1 


-0.04 


NGC4649 


0.65 


±0.12 


-0.05 


NGC6482 


0.52 


+0.19 
-0.23 


-0.07 



+0.03 — 
-0.10 

+0.07 
-0.13 

+0.04 

+0.1 

-0.003 

+0.09 

-0.01 

+0.06 

-0.008 

-0.04 



+0.12 
-0.09 
+0.07 
-0.05 
+0.11 
-0.09 
±0.2 
+0.08 
-0.06 

±0.13 

+0.11 
-0.09 



±0.16 

+0.41 
-0.08 
+0.46 
-0.14 

±0.6 

+0.27 
-0.03 
+0.57 
-0.12 
+0.37 
-0.03 



-0.14 

+0.12 

-0.05 

+0.01 

-0.03 

+0.02 

-0.004 

+0.04 
+0.02 
-0.05 



+0.06 
+0.12 
+0.06 

+0.03 
-0.04 

+0.17 
-0.06 
-0.13 



-0.05 

-0.14 

-0.04 

-0.006 

-0.02 

-0.13 

+0.005 
-0.03 



-0.003 

-0.001 

-0.002 

-0.009 

-0.010 

-0.05 

-0.01 



Note. — The estimated error-budget for each of the galaxies. Excepting the statistical error, these values estimate 
a likely upper bound on the sensitivity of the (best fit) value of each parameter to various data-analysis choices, 
and should not be added in quadrature with the statistical errors. The systematic uncertainties on M v i r and c are 
estimated for the NFW+stars model. In addition to the best-fit values, we show the 90% confidence interval for 
each parameter (Astat). We also show estimated upper-limits on the systcmatics likely to arise by making various 
changes to our default analysis choices. This includes adopting the N04 DM model (AN04), varying the shape of 
the stellar mass component (Astars), varying the background (Abkd), excluding data in the vicinity of asymmetries 
(Aasym), adopting alternative temperature models (Atcmp), changing spectral analysis choices (Aspectra), varying 
the distance (Adist) or assuming that all of the stellar baryons arc in the central galaxy (Abaryons). 



In this section we address the sensitivity of our results 
to various data-analysis choices which were made. An 
estimated upper limit on the sensitivity of our results to 
these choices is shown in Table [| These numbers reflect 
the sensitivity in the best-fit parameter to each potential 
source of systematic error, and we stress they should not 
be added in quadrature with the statistical errors. We 
outline in detail below how each of these systematics were 
estimated. Those readers uninterested in the technical 
details of our analysis may like to proceed directly to 



7.1. DM profile shape 

As discussed above, we experimented with replacing 
the NFW model by the revised N04 model, which is 
less cuspy. We caution that the M v ; r -c relation was de- 
rived assuming NFW. We fixed the a parameter for this 



model to 0.17, the mean value determined from simula- 
tions since the inner slope of the DM halo is degenerate 
to some degree with the stellar mass. The quality of 
the fits (N04+stars, AC N04+stars) were typically sim- 
ilar to those using the simple NFW profile. There were 
some slight systematic differences in the inferred M v j r 
and c as compared to NFW. It is interesting to note that 
this model, which is less cuspy than NFW, gave slightly 
larger M»/Lk, although not sufficiently to bring our mea- 
surements completely into agreement with the synthetic 
estimates. For the adiabatically-compressed N04 model, 
M<,/Lk did not increase, but this is unsurprising since 
the stellar component significantly modifies the shape of 
the inner DM halo in this model. We note that the pre- 
dicted typical inner slope for ACDM halos is still under 
debate. If, instea d of the N04 pro f ile wc had adopted the 
cuspier profile of Diemand et al. ( 2005| ), then the resul- 
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tant M*/Lkwou1c1 have been even smaller, and in worse 
agreement with stellar population models. 

7.2. Shape of the stellar potential 

To account for the stellar component, we adopted an 
H90 model, the effective radius of the model being fixed 
to that determined by 2MASS. However, it is not en- 
tirely clear that the H90 model is an adequate descrip- 
tion of the stellar mass. There are some deviations be- 
tween H90 and the de Vaucouleurs model fitted as the 
de facto standard to the optical light profiles of elliptical 
galaxies, particularly in the critical central regions. Fur- 
thermore, the K-band light profiles of elliptical galaxies 
may, in fact, be better described by the Sersic profile 



jet and AGN in NGC4261. We also excluded data from 
an off-axis X-ray asymmetry in NGC 4 125 and excluded 
data for NGC 4472 outside 6', where Irwin & Sarazin 



(e.g. Brown et al. 2003). To investigate the sensitivity 



of our results to the H90 assumption, therefore, we ex- 
peri mented with adopting a S ersic stellar mass potential 



(e.g. Prugniel fc Simien 1997 ). To determine the two pa- 



rameters of this model (the Sersic index and half-light 
radius) we obtained the K s -band 2MASS images of each 
galaxy from NED, and fitted the surface brightness pro- 
files using dedicated software. A Sersic model fitted the 
K s -band light profile of each galaxy in the radial range 
5"-3' reasonably well. The fitted profiles tended to be 
slightly more centrally peaked than H90, which resulted 
primarily in slightly lower inferred M„«/Lk ratios when 
adopted as mass models. We also experimented with 
r eplacing the H90 mode l with a de Vaucouleurs model 
( Mcllier fc Mathez 1987 ), and adopting R e g values from 
Pahrcj ( |1999| )! 

Elliptical galaxies exhibit radial colour gradients, 
which may reflect gradients in the metallic itv or age of 



the stellar population (see discussion in § |8.4| ). These 
may therefore imply a radial gradient in the stellar M/L 
ratio. It is beyond the scope of this present work to take 
such a gradient into account. However, we investigated 
the sensitivity of our results to the precise shape of the 
optical light profile we adopted by experimenting with 
replacing the K-band R e ff for each galaxy with the (typi- 
cally larger) B-band value listed in RC3. For NGC 6482, 
for which R e g is not listed in RC3, we simply increased 
R cff by 50%. 

7.3. Background subtraction 

One of the major potential sources of systematic un- 
certainty in measuring the mass profiles of galaxies is 
the background subtraction technique. This is especially 
important in the low surface-brightness regime at large 
radii in our galaxies. In order to estimate the likely mag- 
nitude of uncertainty arising from our modellin g, w hen 
initially fitting the background components (§ 3J) we 
artificially adjusted the slope of the instrumental back- 
ground components, which dominate at high energy, to 
the limits of their 90% confidence regions, refitting the 
other components and then refitted all the spectra with 
these revised background models. 

7.4. X-ray asymmetries 

We note that there are some low-level asymmetries in 
the X-ray surface brightness profiles (§ [I]). In order to 
assess the potential impact of these features, we experi- 
mented with excluding or including the features. In par- 
ticular, we tried excluding data from the vicinity of the 



( 1996 ) pointed out that the X-ray data become asymmet- 
ric. These choices did not dramatically affect our results, 
indicating that these features did not indicate a signifi- 
cant violation of hydrostatic equilibrium, at least in an 
azimuthally-averaged sense. To gain an insight into pos- 
sible asymmetries in other sources, we tried re-extracting 
all our spectra, and re-deriving the mass profiles, from 
suitably-oriented semi-annuli (thereby excluding one half 
of the emission from each system). 

7.5. Temperature profile 

In principle multiple temperature profiles may be able 
to fit the same data adequately well but give rise to 
slightly different global halo parameters. In particular 
our constraints upon fb, the computation of which re- 
quires the extrapolation of the density (and hence tem- 
perature) profiles to large radius, may make M v i r and 
c sensitive to this effect. To test this, we experimented 
with cycling through each of our adopted temperature 
profiles (eq Provided the fits were of compara- 

ble quality to our preferred choice, the impact on the 
best-fit parameters reflect the systematic uncertainty in 
this choice. Furthermore, we also experimented with ex- 
cluding the central bin from the temperature profiles of 
NGC 1407 and NGC 4649, which may indicate a central 
disturbance (although there is no obvious X-ray morpho- 
logical disturbance in this region). These choices did not 
strongly affect our results. 



7.6. Spectral-fitting choices 

A variety of choices are made in the spectral-fitting, 
each of which can affect, to some degree, the inferred 
p g and temperature in each radial bin. A thorough dis- 
cussio n of these effects is given in Humphrey & Buote 
p06[ ). 

Column density. In order to take account of possible 
local deviations in the line -of-sight Nji from the value of 
Dickey & Lockman (1990), we experimented with allow- 
ing Vh to vary by ±25%. 

Bandwidth. To estimate the impact of the bandwidth 
on our fits, we experimented with fitting the data in the 
energy ranges 0.7-7.0 keV, 0.5-2.0 keV and 0.4-7.0 keV, 
in addition to our preferred choice 0.5-7.0 keV. 

Plasma code. There are some uncertainties in the cor- 
rect modelling of the individual emission lines, in par- 
ticular those of Fe. This can systematically lead to 
differences in the inferred temperature and density de- 
pending on choice of plasma code. We therefore ex- 
perimented with replacing the APEC model with the 
MEKAL plasma model. 

Unresolved source component. We included a 7.3 keV 
bremsstrahlung component to account for unresolved 
point sources within D25 . This model is generally a good 
fit to the resolved point sources in early- type galaxies, 
but is an empirical result which may not be appropriate 
to model all unresolved sources in all early-type galax- 
ies. We therefore tested the sensitivity of our results to 
this model, by replacing the bremsstrahlung component 
with a simple power law (with T =1.5) or varying the 
temperature of the component by ±25% 
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7.7. Distance uncertainty 

The estimated distance to the object enters into our 
mass determination (Eq [I]) primarily through the impact 
on the radial scale. To assess its impact on our fitting, 
we experimented by varying the distance by ±20%. 

7.8. Stellar baryon fraction 

In our analysis, we restricted ft, to enable interesting 
constraints on M v ; r and c. For the stellar contribution 
to the baryon fraction, we assumed that mass is divided 
among group members following the B-band light, which 
is not formally correct since M*/L ratios are very sen- 
sitive to the age of the stellar population. To estimate 
how much impact this makes to our fits, we experimented 
with assuming that all the stellar mass is in the central 
galaxy, which should place an upper limit on the uncer- 
tainty arising from this choice. 

8. DISCUSSION 

8.1. Hydrostatic equilibrium 

Our fit results provide strong evidence that the gas 
is in hydrostatic equilibrium in these systems. Despite 
highly nontrivial temperature and density profiles, we 
were able to recover smooth mass profiles in remarkably 
good agreement with expectation for these systems, us- 
ing two complementary techniques. If the gas is signifi- 
cantly out of hydrostatic equilibrium, this would repre- 
sent a remarkable "conspiracy" between the density and 
temperature profiles. It is unsurprising that the gas is 
close to hydrostatic equilibrium in these systems, since 
we took care to choose objects with relaxed X-ray mor- 
phology. Based on N-body/ hydrodynamical analysis, 
X-ray measurements are expected to give reliable con- 
straints on the DM in sys tems without obvious substruc- 
ture ( |Buote fc Tsai||l995| ). 



Further support for hydrostatic equilibrium is pro- 
vided by the general agreement between our measured 
M*/Lk ratios and those predicted by SSP models, cou- 
pled with the agreement between the measured M V j r -c 
relation and that expected. Similarly a comparison be- 
tween our results and masses determined from stellar 
dynamics provides even more evidence that our mea- 
sured mass profiles are reliable. Dynamically-determined 
Mgrav/Le within the B-band K c fs are typically fou nd to 
be ~4-10 QGcrhard et alj p00l| ; [Lrujillo ct aL||200^ ). We 
found M, " " 



,_. , R v/Lr with in the B-band R n ff (taken from RC3 
Faber et al. 1989) for our systems ranged from ~3 to 



~8, in good agreement with this result. For NGC4649, 
outside R e ff there is excellent agreement between our 
measured M grav /L profile and that obtained from glob- 
ular cluster kinematics, although at small radii the X- 
ray data lie ^30% lower (K. Gebhardt et al, in prepa- 



ration), van dcr Marel (1991) constructed stellar kine- 
matical models tor 5 galaxies in our sample (NGC 720, 
NGC1407 and NGC4261, NGC4472 and NGC4649), 
under the assumption of a constant M/L profile. Strictly 
speaking a direct comparison cannot be made between 
their M grav /LB measurements and our results since our 
data indicate this assumption is incorrect. However, if we 
simply assume that these M/L ratios represent those in- 
tegrated out to R e ff, the X-ray inferred masses vary from 
~40 % lower to ^10% h i gher than those from kinemat- 
ics. Kronawitter et al. (2000) report M grav /LR~8±1.5 



for NGC 4472 within ^50", at which radius our X-ray 
determined value is ~50% lower. The discrepancies be- 
tween the X-ray and dynamical masses are only modest 
(the X-ray mass being on average ~20% lower), indicat- 
ing that the data must be close to hydrostatic equilib- 
rium. Turbulence is expected to contribute only ~10% 
pressure support in clusters, which are believ ed to be 
more turbulent than galaxies, ( Rasia et al. 2006 ). There- 
fore, on a case-by-case basis, the observed differences are 
most likely a manifestation of the mass-anisotropy de- 
generacy (e.g. Dekcl et al. 2005 ). 



8.2. Mass profiles 

We obtained detailed mass profiles for 3 galaxies and 
4 group-scale systems, out to ~10R c ff. The data clearly 
show M/L profiles which are ^flat within R c ff and rise 
considerably outside this range. This confirms the pres- 
ence of substantial DM in at least some early-type galax- 
ies and indicates that a stellar mass component domi- 
nates within ^R c ff. This is consistent with studies of 
stellar kine matics and similar to the ma ss decomposition 



analysis of Brighcnti fc Mathews ( 1997 ) 



The data are well-fitted by a model comprising a stel- 
lar mass (H90) component and an NFW DM profile. 
Omitting the stellar mass component led to systemat- 
ically poorer fits, smaller M vir and vastly larger c (^20), 
which are inconsistent with the predictions of ACDM. 
This effect is easy to understand — if we add a compact 
stellar mass component to an (extended) NFW profile, 
we increase the mass in the core which, by definition, 
makes the halo more concentrated. However, i t is not 
entirely c l ear w hether this effect, pointed out by Mamon 
Lokag ( 2005 ), can completely acco unt for the signifi - 
cantly steeper M v i r -c relation found by 5ato et al. ( 2000]). 
Based on o ur analysis of group-scale halos jOastaldcll~ 



it al.| 2006 ) we found that the inclusion of the stellar 
mass component does not have a strong effect on c in 
most systems with M v ; r > 2xl0 13 M Q , provided the data 
are fitted to a sufficiently large fraction of R v i r - The 
data did not allow us to distinguish statistically between 
the simple NFW+stars model and scenarios in which the 
DM halo experiences adiabatic compression due to star 
formation (however, see § S.4), or the NFW profile was 



replaced with the alternative N04 profile. 

Comparing our inferred M v ; r and c to the predictions 
of ACDM we find general agreement. There is some ev- 
idence, however, that the concentrations are systemati- 
cally higher than one would expect, although the error- 
bars are typically l arge. Such a trend is als o seen in our 
analysis of groups (Gastaldcllo ct al. |2006 ). Whilst the 
slope of the M V j r -c relation therefore implied by our data 
is difficult to explain by varying the cosmological pa ram- 



eters within reasonable limits (Buote et al. |2006a ) , we 
suspect that the discrepancy can be resolved by taking 
into account the selection function of our galaxies. 

Our data were not selected in a statistically complete 
manner and, by choosing objects with relatively relaxed 
X-ray morphologies we are probably selecting objects 
which have not had a recent major merger. This sys- 
tematically biases us towards early-forming, hence higher 
concentration halos. In fact, it is striking that all three 
de facto ga laxi es in our sample are relatively isolated 
systems (§ |S.3| ). Such systems pr eferentially might b e 
expected to occupy high-c halos ( ^entner et al. 200H ), 
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which does appear to be the case for 2 out of 3 of the 
galaxies. We will return to these issues in detail in Buotc 
et al| (|2006a|) . 

8.3. Galaxies, Groups and Fossil Groups 

All three of the lowest-mass systems in our sample 
are very isolated optically. NGC 6482 matches the isola- 
t ion criteria adopted to id entify so-called "fossil groups" 
flKhosroshahi et al.| p004| ) . NGC 4125 and NGC 720 are 
both listed as "groups" in G93, but closer inspection 
actually reveals they are also very isolated. Excepting 
the central galaxy, only one of the putative members of 
the NGC 720 "group" listed in the G93 catalogu e (which 
omits the dwarf galaxy population studied by Dressier 
t al. 1986), actually lies within the projected R V i r (but 
outside 0.75xR v j r ) and it is 2.4 magnitu des fainter in B 



E 



than the central galaxy. Dressier et al. remarked upon 



the optical isolation of this galaxy. Of the two putative 
companion galaxies to NGC 4125 given in G93 which lie 
within the projected R V j r (but outside 0.67xR v i r ), both 
are much fainter (by 2.3 and 3.9 magnitudes, respec- 
tively) in B than the central galaxy. In contrast, the 
four remaining s ystems in our sample a re much less op- 
tically isolated. Schindler et al. (199E) show the clear 
over-density of early-type galaxies around NGC 4649 and 
NGC 4472, and almost 60 gro up members are associated 
with these systems by G93. Gouldj (1992) identified at 
least 10 members of the NGC 1407 group, from the dy- 
namics of which he inferred a m as s broadly consistent 



with o ur measured M V i r (§ 8/7). Helsdon fc Ponman 



( [2003D 



report 57 galaxies associated with the NGC 4261 
group within ~1 Mpc projected radius, which is consis- 
tent with our measured R v i r . 



Rat her than an isolated galaxy Khosroshahi et al 
( p004D identify NGC 6482 as a "fossil group". Fossil 
groups are group-sized X- ray halos centred on essentially 
a sing l e elliptical galaxy ( Ponman et al. 1994 ; Vikhlinin 
|st al. 1999; Jones et al. |2005 ). The typical interpre- 
tation of these objects is groups in which all of the 
L* members have merged. Confusingly, using almos t 
the same selection criteria, O'Sullivan fc Ponman| ( 2004| ) 



classify the galaxy NGC 4555 as an "isolated elliptical 
galaxy" and posit a very different formation scenario. 
This object appears to be more massive than NGC 6482; 
the authors found M grav ~ 3 x 1O 12 M0 within 60 kpc 
which, assuming an NFW profile with c=15 would imply 
M v i r ~ 2 x 10 13 M Q . Nonetheless, both of these systems 
have more in common (both optically and in the X-ray) 
with each other, and the other isolated ellipticals in our 
sample, than, for example, the massive (M v ; r > 10 14 M^). 



hotte r (kT^ 2 keV) fossil groups considered by |Vikhlinin 



|st al 



( 1999 ). We suspect that the distinction made be- 
tween "isolated elliptical" and "fossil group" for these 
two systems is largely semantic, and consider NGC 6482 
more properly an isolated galaxy, too. 

The clear division in the galaxy content of our sample 
clearly lends itself to the nomenclature "galaxies" for the 
three lowest-mass systems, and "groups" for the others. 
Strikingly, this separation between galaxies and groups 
also appea rs c onsistent with a difference in temperature 



profiles ( 
surate with ]V 



(§M) 



That this distinction appears commen- 
r ~ 1O 13 M0 is suggestive that this mass- 
scale may be a useful yard-stick with which to compare 
to other systems. The error-bars on our mass estimates 



are sufficiently large that the 90% confidence regions of 
several of the objects (notably NGC 720, NGC 6482 and 
NGC 1407) actually straddle 10 13 M Q . However, it is 
clear that on average, the systems with M v i r < 10 13 M Q 
are galaxies. We note that the M v i r adopted here is that 
before any tidal truncation which is almost certainly oc- 
curring as NGC 4472 and NGC 4649 merge with Virgo 
(their untruncated R v i r would stretch much of the dis- 
tance to M87). M v ; r does not exactly correlate with for- 
mation epoch, so that lower-mass halos may still be in 
the process of forming (hence contain multiple galaxies of 
similar magnitude) , and more massive halos may contain 
single, dominant ellipticals (fossil groups). Nonetheless, 
classifying halos primarily on the basis of M v ; r provides 
a straightforward way to locate them in the formation 
hierarchy. Traditionally, galaxy-like and group-like sys- 
tems are distinguished on the basis of local over-densities 
of galaxies. However, placing optically-identified groups 
into a cosmological context requires a firm understanding 
not only of the formation of DM halos but also how galax- 
ies p opulate them, which is much less well- understood 
(e.g. Kravtsov et al. 2004 ). This problem is compounded 
by the difficulties faced by optical group- finding algo- 



rithm s in identifying very poor groups (e.g. Gerke et al 
2005). Not only can a significant fraction of putative 



groups be chance superpositions of galaxies, particularly 
along filaments, but adjacent groups can be merged, such 
as happened for NGC 4649 and NGC 4472 in G93. If 
there are only a few identified members, small-number 
statistics and the tr eatment of in terlopers can affect their 
interpretation (e.g. Gould 1993j ). To complicate matters 
further, some authors refer to any over-density of galax- 
ies as a group, even a Mi lky W ay-sized galaxy and its 
dwarf satellites (e.g. |Tullyj ^005|) . 



8.4. Stellar Mass-to-Light Ratios 

Comparing our measured stellar M/L ratios to the 
predictions of simple stellar population (SSP) models, 
w e found r eason able agreement provided one assumes 
a Kroupa ( |2001 ) IMF. There is modest disagreement, 
even when the less-cuspy N04 DM model was adopted. 
Considering the uncertainties in the SSP modelling (dis- 
cussed below), however, we believe this discrepancy is 
not significant. If we allowed the DM profile to be modi- 
fied by adiabatic compression, we obtained substantially 
smaller M»/Lk values from our data, (since it increases 
the cuspiness of the halo) which are more discrepant with 
the SSP models. This result casts doubt on AC being as 
significant an effect as currently modelled. However, the 
data alone did not allow us statistically to distinguish 
between the NFW+stars and AC NFW+stars models. 
Nonetheless, this result is joining a growing body of lit- 
erature which similarly calls into question whe t her AC 



opera tes as predicted ( Zappacosta et al 
' t al.||2006| ; |Sand et al.|[2004| ) 



2006j ; [Kassin 



There are a number of major uncertainties in the com- 
putation of the stellar mass-to-light ratios from the SSP 
models. Specifically, the results are very dependent upon 
the assumed IMF, which is not confidently known in 
early- type galaxies. Furthermore there is some evidence 
that early-type galaxies frequently contain multiple stel- 
lar populations of diff erent ages, including a significant 



youn g population (e.g. Rembold et al. 2005 ; Nolan et al 
2006). Depending on the mass fraction of the young 
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component, this may substantially reduce M*/Lk in the 
galaxy, hence possibly reconciling the data and the AC 
NFW+stars model. A small amount of star formation 
may also give rise to a population of stars which can 



can dominate over radiative cooling to produce a nega- 
tive temperature slope. In contrast, if hotter (~l-2 keV) 
baryons are allowed to flow in, radiative cooling domi 



dominate the light in the galaxy core, giving rise to sig- 



nificantlv lower synthetic M>/Lk t han measured. This 
may be the case in NGC 720 (see § |3lj| ). More problem- 
atically, there are known to be significant abundance, 
or possible age, gradie nts in the stellar populations of 



nates to produce a po sitive temperature gradient (Math- 
ews fc Brighcntj|2005j ). It is by no means clear, however 



earl y-type gal a xies (e.g. Tragcr et al 
Arin joto|[l999| ; |Rembold et all ^005|) 



2000| [Kobayashi & 
which would trans- 
late into stellar M*/Lk gradients. Our simple modelling 
did not allow us to account for such an effect per se. 
Although we suspect that such gradients will primarily 
lead to a M*/Lk value which reflects an average for the 
galaxy, M„/Lk does depend to some extent upon the 
shape of the assumed stellar potential. Properly taking 
account of this effect is beyond the scope of this present 
work, but may bring the synthetic M/L ratios and our 
results into better agreement. Clearly this is only one 
of a number of other systematic effects which may also 
reconcile the slight discrepancy (Table |). 

8.5. Baryon fractions 

An interesting result from our analysis is that these 
systems, despite having masses > 5xlO 12 M , do not ap- 
pear in general to be baryonically closed. To some ex- 
tent this trend was enforced by applying Eq |^ to con- 
strain the data. However, the excellent fits we obtained 
by this method, in conjunction with the good agreement 
between the measured M v ; r -c relation and the predictions 
of ACDM and, crucially, our measurements at the group 
scale (whic h do not employ this restriction: Gastaldcllo 
[ t a,L||2006| ), indicate that the inferred f b (~0.04-0.09; Ta- 
ble"!) are accura te. Furthermore, if we relaxed this con- 
straint and instead restricted fb to a finite range, we also 
found that the data tended to favour modest values of fb- 
In particular, for any given system, the measured M v j r 
and fb were strongly anti-correlated, so that our upper 
M v i r constraint is in part imposed by the lower limit we 
place on fb. Given the shapes of the M v i r -c contours 
(Fig , it is clear that good agreement with the M v i r -c 
relation predicted from simulations tends, therefore, to 
require rather modest values of fb- This would suggest 
that strong feedback plays an important role in the evo- 
lution of these objects. 

8.6. Temperature profiles 

By inspection of the temperature profiles (Fig 0) it is 
immediately clear that, for all of the galaxy-scale sys- 
tems in our sample the temperature profiles have nega- 
tive gradients. In contrast the group-scale objects have 
positive temperature gradients, similar to observations 
of other X - ray bright groups and clusters ( Gastaldcllo 
|st al.| [2006| ; |Vikhlinin etaT] |2005| ; |Piffaretti et al.| p005| ) . 
This radical difference in the temperature profiles seems 
consistent with our division of galaxies and groups at 
M v i r ~ 10 13 M Q . The origin of this distinct demarca- 
tion between objects around 10 13 Mq is unclear, however. 
Negative temperature gradients are expected for isolated 
galaxies containing relatively cool gas, such as that aris- 
ing from stellar mass-loss. In the deep stellar potential 
well, compressive heating of the gradually inflowing gas 



why the hot baryons appear to be present only in the 
systems with M v i r > 10 13 M Q . One possibility is the local 
environment; all of the galaxy-scale objects are rather 
isolated, whereas the groups NGC 4472 and NGC 4649, 
in particular, are found in a relatively dense cluster envi- 
ronment, which could provide a reservoir of hot baryons. 
However, such an explanation cannot easily account for 
the positive temperature gradient in NGC 1407, which is 
comparatively isolated, or the isolated system NGC 4555, 
which appears only slightly more massive than our galax- 
ies. 

It is possible that selection effects may have played 
some role in the bimodal temperature profile behaviour, 
since both NGC 4125 and NGC 6482 are classified in 
NED as LINERS, and NGC 720 has a dominant young 
stellar population (Appendix [A]) . However, none of these 
systems show strong X-ray morphology disturbances in 
the core, which might indicate a substantial energy input 
from star-formation or AGN activity. In any case, the 
cooling time in the core of NGC 720 is only ^200 Myr, 
substantially less than the implied time since the last 
major burst of star formation, and so the negative tem- 
perature gradient cannot simply be related to energy in- 
jection during a starburst. Furthermore, at least two 
of the group-scale systems also harbour AGN and do not 
show obvious negative temperature gradients in the core. 

Another example of an object we believe to be a galaxy 
(rather than a group) which exhibits a negative tem- 



erature gradient is the SO NGC 1332 ( Humphrey et al 
004 ) . A possible counter-example to this trend might by 



the "isolated elliptical galaxy" NGC 4555, which exhibits 
a temperature profile akin t o the groups in our sample 
( O 'Sullivan fc Ponman 2004). However, as we discuss in 
8.3, this probably has comparable M V i r to the groups. 



Another intriguing feature of two of the group scale 
objects is a central temperature p eak, similar to a fea- 
ture we found in the cluster A 644 (Buote et al. 2005). In 
that system, we found a significant offset between the X- 
ray centroid and the emission peak in an otherwise fairly 
relaxed object. We suggested that both of these features 
may be related to the cD "sloshing" in the potential well 
of the cluster, which is relaxing following disturbance by, 
for example, a merger. We do not find obvious evidence 
of a similar offset in either NGC 1407 or NGC 4649. How- 
ever, these groups may be in a comparably more relaxed 
(evolved) state than A 644. Alternatively, the central 
peaks may be related to past AGN activity heating the 
gas in the core of the galaxies, from which the system has 
had time to relax dynamically but not cool completely. 

8.7. Is NGC1407 a "dark group"? 



Bas ed on the group member dispersion velocity Gould 



(|1993|) suggested that NGC 1407 may lie in a massive 
(>a few xl0 13 M Q ) DM halo. Although such a conclu- 
sion was strongly dependent on the association of the 
galaxy NGC 1400, which exhibits a large peculiar veloc- 
ity, with the group, we can now confirm the presence 
of a substantial DM halo around this system. Both the 
temperature profile and our best-fit mass are similar to 
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the bright X-ray group NGC 5044 flBuote et al] j2006b|) 
and yet it is almost 2 orders of magnitude fainter in Ly 



NGC 5044 appe ars to be close to baryonic closure ( Math- 
ews k t al, 2005 ), and so has likely retained most of its 
On the other hand NGC 1407 is not 
0.06) and so the loss 



large gaseous halo, 
baryonically closed (we estimate 
of much of its hot gas envelope easily explains its lower 
Lx/Lb- Since the masses of the two systems are not 
considerably different, this points to substantial varia- 
tion in the evolutionary history of these two groups. In 
particular, feedback may have operated more efficiently 
in evacu ating the gas from NGC 1407. 

Gould's preferred mass estimate (~ 10 14 M Q ) would 
imply a remarkably high M/L ratio for the system 
(M vir /L B ~9OOM /L ), making NGC 1407 a bona fide 
"dark group" . The existence of such an object would 
provide a valuable insight into the process of star forma- 
tion in DM halos, as it would imply star formation was 
somehow inhibited in that system. This mass estimate 
is, however, considerably larger than our preferred value 
~ 1.5 x 10 13 M Q , which implies a more modest M/L ra- 
tio (M V i r /LB~14OM0/L Q ). To some extent, though, our 
constraint on fb, which was necessary to obtain inter- 
esting M v ; r constraints, has probably enforced this be- 
haviour. Such a restriction may not be valid in a sys- 
tem with an unusual star-formation history and so we 
experimented with freeing fb. To enable M v ; r to be con- 
strained, we rest ricted c to lie o n the best-fit M v i r -c re- 



lation found by Bullock et al. ( 2001 ). The best-fitting 



mass, M vir = {9-7 t l 7 2 S ) x 10 13 M Q , was in good agree 



ment with Gould ( |1993 )'s values, but implies a baryon 
fraction only of ^0.003. Since this fit was statistically 
indistinguishable from the preferred model, we cannot 
determine which mass estimate is more likely. 

9. SUMMARY 

Using Chandra we have obtained detailed mass profiles 
centred on 7 elliptical galaxies, of which 3 were found to 
have de facto galaxy-scale halos, with M v ; r < 10 13 M Q , 
and 4 had group- scale (10 13 M Q <M vir < 10 14 M Q ) halos. 
These represent the best available data for nearby objects 
with comparable Lx- In summary: 

1. The M/L ratio profiles were ~flat within R g and 
rose sharply outside this region, indicating substan- 
tial DM in all 7 systems. 

2. The data were well-described by a two compo- 
nent model, comprising an NFW potential for the 
DM and a H90 stellar mass model. We were not 
able statistically to distinguish between this sce- 
nario and one in which the DM profile was mod- 
ified by "adiabatic compression" due to baryonic 
infall. Similarly, we could not distinguish between 
the NFW and the revised N04 DM halo profiles. 

3. The distribution of the galaxies in the M v i r -c plane 
was in broad agreement with the predictions of 
ACDM, although with a slight trend toward more 
concentrated halos, in good agreement with our 
modelling of X-ray brigh t groups and poor clusters 
(Gastaldello et al, 2006[ ). This probably represents 
a galaxy selection bias to earlier- forming systems, 
and we will d iscuss how we might account for it in 
Buotc et al. (2006a). Allowing AC to modify the 



shape of the DM halo did not appreciably affect 
the M V i r -c relation. 

4. Omitting the stellar mass component resulted in 
systematically poorer fits, smaller M v ; r and un- 
physically large c , confi rming the conclusions of 



Mamon & Lokat (2005). This may explain very 



large values of cfound by som e previous X-ray ob- 



2004). 



serve rs (e.g. |Sato et ajj pOOOt |Khosroshahi et al 



5. For the NFW+stars model, M*/Lk was found to be 
in approximate agreement with the predictions of 
simple stellar popula tion synthesis models, assum- 
ing a |Kroupa| (|200l| ) IMF. The AC NFW+stars 
models have significantly lower M»/Lk which 
seems to cast doubt on the AC scenario, although 
this conclusion is sensitive to the considerable un- 
certainties in the theoretical modelling. 

6. Despite having M vir > 5x 10 12 M Q , typically 
fb~0.04-0.09 for each galaxy, implying that feed- 
back has played an important role in the evolution 
of these systems. 

7. The temperature profiles of the galaxy-scale 
systems all exhibited negative radial gradients, 
whereas the group-scale objects exhibited positive 
gradients, similar to the "Universal" temperature 
profiles being found in other X-ray bright groups 
and clusters. This implies a strict line of demarca- 
tion between systems at M v i r ~ 10 13 M Q . 

8. In two of the groups, we found central temperature 
peaks, similar to t hat found in the cluster A 644 
( Buote et al. 2005 ), but no obvious central distur- 
bances in X-ray morphology. This may relate to 
past AGN activity, following which the heated gas 
in the core of the galaxy has relaxed but not cooled. 



9. We confirm the suggestion of Gouldj ( 1995 ) that 
the elliptical galaxy NGC 1407 lies at the cen- 
tre of a massive DM halo, possibly making it a 
"dark group" with an unusually large M/L. Our 
best-fi tting; M v i r is considerably lower than that of 



Gould, implying M/L more consistent with normal 



groups. Nonetheless, if we relax the assumptions of 
our modelling very large masses (M v ; r ^ 10 14 M Q ) 
are allowed. 
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TABLE A7 
Stellar population parameters 



Galaxy 


indices 


rcf. 


[a/Fe] 


age 


[Z/H] 


< [Z/H] > 










(Gyr) 






NGC 720 1 " 


H/3, Mgb, Fe5270, Fe5335 


2 


0.37±0.05 




0.65±0.13 


0.48±0.18 


NGC 1407+ 


H/3, Mgb, Fe5270, Fe5335 


1 


0.33±0.02 


12±2 


0.35±0.06 


0.08±0.06t 


NGC 4125 


H/3, Mgb, Fe5270, Fe5335 


3 


0.33±0.16 


13±8 


0.16±0.25 


-0.11±0.25} 


NGC 4261 


H/3, Mgb, Fe5270, Fe5335 


2 


0.25±0.02 


15±1 


0.30±0.03 


-0.03±0.10 


NGC 4472t 


H/3, Mgb, Fe5270, Fe5335 


2 


0.25±0.03 


9±2 


0.36±0.05 


0.17±0.12 


NGC 4649 


H/3, Mgb, Fe5335 


2 


0.25±0.02 


13±2 


0.41±0.04 


0.23±0.15 


NGC 6482 


Mgb, Fe5270, Fe5335 


3 


0.30±0.15 


12 


0.28±0.15 


0.06±0.15t 



Note. — Stellar population parameters determined from Lick index fitting. The indices used 
in fitting are listed (indices), as is the reference (ref) whence they were taken. Those mean stellar 
abundances (< [Z/H] >) marked | were estimated from the central ahundanro ([7, /H]p) adopting 

the mean abundance gradient [Z/H]o — < [Z/H] > = 0.27(sce Humphrey ^ Bugle ZDQ£j ) ^Vhere 

the parameter was frozen . Table references 



no error-bar 


is eiven 


Tra.gnr et al 


1 1' 


Humphrey & 


Buotc] 



3= 



(20oe) 



Trager et al 



(1998) Results for galaxies marked 



ticuing et 



aT 



I gUOgj ), 2— 
were taken from 



under grant NNG04GE76G issued through the Office of Space Sciences Long-Term Space Astrophysics Program. 

APPENDIX 

STELLAR POPULATION PARAMETERS 

The mass-to-light ratio of a stellar population is dependent upon both the age and the metal abundance ([Z/H]) of 
the stars. To estimate these quantities we searched the litera ture to obtain Lick/I DS indices for each galaxy, which 



we fitted with the s i mple stellar population (SSP) models of Thomas et al. (2003), using the technique outlined in 
Humphrey & Buotc (2006). Briefly, we constructed a model by linearly interpolating the SSP models as a function 
of stellar age, metallici ty and q-elemcnt to Fc ratio, which was then fitted via a \ 2 minimization technique to those 
indices shown in Table A7. Trager et al. ( [2000 ) pro vided indices measur e d in t wo apertures, which enabled us to take 
account of any abundance gradients, as outlined in Humphrey fc Buotc] (2006). Where only a central Lick index was 
available, we estimated the total emission-weighted abundance by correcting the central metallicity by -0.27 dex. We 
did not attempt to take account of possible age gradients. T he re sults, including the Lick indices adopted and the 
reference whence the indices were obtained, are shown in Table A7. This method implicitly assumes that all the stars 
were created in a single burst of star form ation, which may be o ver-simplistic if th ere are, in fact, multiple bursts of 



star formation in early-type galaxies (e.g. Rcmbold et al. 2005 ; Nolan et al. 2006] ). We note that we were not able 
to obtain acceptable solutions for NGC 6482 if we used the H/3 index, which is the most sensitive age indicator. This 
galaxy is classified in NED as a LINER and is rather blue for an old stellar pop ulation (Mb — M k = 3.4 ± 0.2, whereas 
a 12 Gyr, solar abundance population is expected to have Mb — Mk — 3.9: Maras ton] I1998T). Both of these factors 



might suggest the presence of a significant young population of stars (although sec Cid Eernandes et al. 2004 ). It is 



beyond the scope of this present work, however, to attempt to take account of this effect. 
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